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ABSTRACT 


Measurements of the longitudinal sound velocity Vi 
in hep He* Single crystals were made as a function of direction 
relative to the hexagonal axis for molar volumes between 20.3 
and 17.3 cm?/mole. The direction was determined by measuring 
the birefringence of the crystals along the sound path. The 
variation of Vi with molar volume is much stronger than 
predicted by theory but agrees with neutron scattering data. 
Agreement is also found with theoretical predictions for the 
anisotropy of the velocity. 

All five elastic constants were determined from these 
measurements and from published compressibility data. They 
are compared with existing theoretical and experimental 
constants for solid helium. 

The temperature dependence of Vv, was measured at a 


molar volume of 20.42 cm/mole and was related to the tempera- 


ture dependence of the adiabatic compressibility. 
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1. INTRODUCTION 


Solid helium attracts the interest of experimental 
and theoretical physicist for various reasons. It can be 
prepared. with a chemical and isotopic purity that is 
impossible in other substances. Its density can be doubled 
with modest equipment and finally some hope still persists 
that the spectacular behaviour of liguid helium below a 
temperature of 2 K would be reflected in some properties 
of the solid. From a,theoretical point.of view, the 
challenge consists in the complete breakdown of classical 
lattice dynamics when it is applied to solid helium, most 
drastically demonstrated by the imaginary sound velocities 
that classical theory predicts. 

The activity shown by theorists in the past ten 
years to remedy the situation has unfortunately not been 
matched by many experiments on sound propagation. | Measurement 
of the longitudinal sound velocity in solid Ee and oie by 
Vignos and H. Fairbank (1961, 1966) on polycrystals gave 
values with a spread of 12%, indicating the presence of 
large cristallites or possibly single crystals and anisotropies 


of that order. The experiment covered a pressure range of 
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26 to 140 atm, but nothing definite could be concluded about 
the volume dependence of the velocity because the orientation 
of the crystals was unknown. Shortly afterwards, Lipshultz 
and Lee (1965) observed transverse sound in He* between 26 
and 29 atm and found a somewhat larger anisotrophy in the 
hexagonal phase. Longitudinal sound in bcc He? was also 

seen by Abel et al. (1961). 

From the theoretical side, Nosanow and Werthamer 
(1965) found an anisotropy in the longitudinal velocity that 
was. larger than the observed 12% and also a smaller volume 
dependence. At that time it was thought that an improved 
theory (Gillis et al, (1968)) would be accurate to within a 
few percent. A meaningful comparison with experiment requires 
therefore that the sound velocity is measured as a function 
of orientation so that it can be clearly decided whether any 
discrepancies are real or simply due to anisotropy. 

Another incentive was provided by a theoretical paper 
by Guyer. (1966) on sound absorption in solid helium, which 
led to the observation of second sound by Ackerman et al. 
(1966). It predicted two peaks in the sound absorption as 
a function of temperature, one associated with normal pro- 
cesses.,. the other with umklapp processes. To observe sound 


attenuation turned out to be difficult and no reproducible 
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results have been obtained so far (apart from the rule that 
the absorption increases with decreasing temperature). We 
found it therefore advisable to first gain more experience 

in ultrasonics and in preparing single crystals of solid 
helium. The present work should therefore be regarded as a 
first step in a project to measure sound absorption. Shortly 
after starting it, the birefringence of hexagonal close packed 
(hep) He* was measured by Vos et al. (1967) and by Heybey 

and Lee (1967). This offered a convenient method of finding 
the orientation of helium crystals. At the same time I 
realized that the sound velocity in hexagonal crystals depends 
only on the angle between the crystallographic c-axis and the 
propagation direction. Therefore, one does not have to know 
the direction of all the crystallographic axes but only the 
c-direction relative to the direction of sound propagation, 
and this is exactly what the birefringence technique pro- 
vides. 

The impact of these facts led me to abandon the 
absorption measurements and to measure the longitudinal sound 
velocity as.a function of orientation and later also as,a 
function of pressure and temperature. Simultaneously but 
independently, the orientation dependence of both longitudinal 


and transverse velocity was measured by a group working under 
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Professor D. M. Lee at Cornell University. They grew the 
crystals from superfluid helium and also use the birefringence 
technique to find the crystal orientation (Heybey and Lee 
(1968), Crepeau et al. (1968), Lee et al. (1969)). 

Preliminary results of our work have already been 
reported elsewhere (Franck (1968), Franck and Wanner (1969), 


Wanner and Franck (1970)). 
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2. THEORY 


Much theoretical work has been done both on sound 
propagation, lattice dynamics and solid helium, so much. in 
fact that even of the review papers only a small portion will 
be mentioned here. 

Classical lattice dynamics is covered extensively in 
Born and Huang's book (1962) and by Maradudin, Montroll and 
Weiss (1963), but I. have found Leibfried's handbook article 
(1955) the clearest introduction to the subject. Anharmonic 
effects are treated by Leibfried and Ludwig (1961). About 
30 papers on lattice dynamics of solid helium have been 
published since 1962 and are reviewed by Werthamer (1969). 

I have made heavy use of this paper in the paragraph on solid 
helium theory. The latest review of experimental data is in 
Wilks book (1967) and in Keller's book (1969). An intro- 
duction to elastic constants. that I have often used is the 
article by Huntington (1958) and the books of the series 


"Physical Acoustics" edited by Mason. 
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2.1 Classical lattice dynamics 


The problem in lattice dynamics is to calculate the 
macroscopic elastic and thermal behaviour of a, solid from the 
given microscopic properties of the particles from which it is 
formed. We will see later that this problem cannot be solved 
in the case of solid helium using classical mechanics alone, 
but since some steps are identical in a classical anda 
quantum mechanical treatment I will give a short outline of 
classical lattice dynamics, based on Ziman (1965). 

To find the equations of motion of each atom we write 


down the kinetic energy of the whole crystal 


where m is the mass of one atom and u the displacement of 
1 
the atom from its equilibrium position 1. The potential 


energy of the crystal is expanded in terms of the cartesian 


components uy of the displacements 


If the approximation is terminated after the second term it 


is called a harmonic approximation. The equations of motion 
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are found using the Lagrangian formulas 


ai if 
mii; =- fF ie aes as 993 
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The: right hand terms are the forces acting on atom 1 due to 
the displacement of atom 1° and can be evaluated if the 
potential between these atoms is given as was assumed in the 
beginning. They clearly depend only on the relative position 
h = 1- 1° of these atoms. 

The periodicity of the lattice reduces the number of 
equations and restricts the solutions to waves because Bloch's 


theorem relates the motion of the atoms at different sites 
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q is the wave vector of a particular solution of the equations 
of motion. The equation of motion for that particular solution 


and a particular atom is then 


mu(q) = (Z2G(h)e"4-. u(gq) = G(g)u(g) 2.7 
G(q) is called the dynamical matrix and is simply the Fourier 
transform of the force constand tensor G(h). To solve that 


equation, we try a solution of the form 


plug it into 2.7 and get 3 simultaneous equations for the 


components ud (depending on the wave vector q) 


The secular equation (that is the condition that the 
determinant of the system of equations is zero) will give 3 
solutions for w(q), corresponding to the three different 
polarizations of the waves. 

In principle then, this is the recipe for calculating 
sound velocities from interparticle potentials: First write 


down the tensor G(h), involving the second derivative of the 
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interparticle potentials of as. many neighbours as-is necessary. 
Next calculate the Fourier transform G(q) for that wave-vector 
g for which the sound velocity is desired. The secular 
determinant of that matrix gives a cubic equation for the 
frequency w. The sound velocity for the wave in the direction 
quis v,= w/,/q]. 

It is customary to restrict the use of the word 
"sound" to the limiting case g>0, that is to long wavelengths 
(long compared with the lattice spacing) and to speak of 
phonons for general g. The elements of the dynamical matrix 
in the long wave length limit are intimately related to the 


elastic constants c. 


she (see equation 2.28). 


2.2 Lattice dynamics of solid helium 


The classical theory of lattice dynamics fails 


dramatically for solid helium. Using a Lennard-Jones potential 
vir) = 4e ((c/r)t? - (o/r)®) 2.10 


for the attraction between two He atoms with parameters e and 
o derived from measured virial coefficients of helium gas it 
predicts a molar volume of 10 em? which is too small by a 


factor 2, a compressibility too small by a factor 30 and sound 
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velocities too high by a factor 4 if compared to experimental 
values at T = OK and at p = 26 bar, the minimum pressure re- 
quired to solidify 9 

An attempt by de Wette and Nijboer (1965) to calculate 
the phonon frequencies classically yielded imaginary frequen- 
cies and sound velocities for volumes. larger than 12 cm?/mole, 
although the lattice is stable for volumes up to 21 cm?/mole. 

The reason for this unusual behaviour is that the 
helium atoms have a large kinetic energy from the zero-point 
motion, an energy comparable or larger than the potential 
energy which has only a very shallow attractive well due to 
the tightly bound electronic shell. This zero-point energy 
can be estimated using the uncertainty principle. Confining 
an atom to its lattice site will define its position with an 
uncertainty of the order of the lattice constant a which will 


lead to an uncertainty in momentum and thereby to a kinetic 


energy per atom of 


a) (Ap)? v h* : Pei eB 


Because of that motion the atom will no longer be at the 
minimum of the potential well and the lattice will expand to 
such an extent that helium solidifies only under pressure, 


even. at 0 K. 
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Using these arguments F. London (1954) has calculated 
the molar volume of the solid at 0 K and found agreement to 
within 15% of the measured value and his simple theory repro- 
duces the volume dependence of the groundstate energy or the 
compressibility remarkably well. 

A more detailed theory of solid helium faces a number 
of problems. Because of the large vibration amplitudes the 
forces on an atom will be highly anharmonic. Figure 1 shows 
a He atom in the potential of its two nearest neighbours in 
a Simple, one dimensional model. Instead of being ata 
minimum, the central atom sits at a maximum and assumed vi- 
brations around this position will be unstable since the second 
derivative of the potential is negative and the vibration 
frequencies therefore imaginary. This is the explanation for 
the imaginary phonon frequencies obtained by de Wette and 
Nijboer (1965). The problem is therefore that the anharmonic 
terms are so large that the quasi harmonic approximation does 
not exist any longer. 

How to introduce into the theory the correlation be- 
tween the movement of neighbouring atoms, the socalled short- 
range or dynamic correlations poses another question: One 
might expect that the large amplitude vibration of. the atom 


would lead to many collisions between neighbouring atoms, 
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involving the strong hardcore repulsion, but in fact the atoms 
will correlate their movements to avoid such collisions. 
Nosanow (1966) points out that correlation lowers both the 
kinetic energy by increasing the radius of the effective 
potential well and the potential energy by increasing the 
probability that two particles are at a distance Cin COrren 
sponding to the minimum in the Lennard-Jones potential. 

Finally, if the theory is formulated in terms of 
phonons it is not obvious a priori whether the large anharmon- 
icity will not lead to unreasonably short phonon lifetimes, 
that is whether phonons do exist in solid helium or not 
(Gillis and Werthamer (1968)). 

Two different approaches have been tried to deal with 
these problems. One is formulated in terms of single particle 
wave functions while the other describes the crystal in terms 
of phonons. 

The basic idea of the single particle picture (Nosanow 
(1964), Fredkin and Werthamer (1965))is to do a time dependent 
Hartree calculation assuming a trial wave function of the so 


called Jastrow form (Jastrow (1955)) 
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FIGURE l 


A helium atom in the Lennard-Jones potential V(r) of its 
nearest neighbours. 


V(r) = 4e((o/r)?* - (c/r)®) 
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The i are single particle wave functions localized around 

the lattice sites 1; while Bive-up) describes the correlations 
between particle j and k. The correlation function f(r) should 
go to zero for r smaller than the hardcore radius, thereby 
decreasing the probability that two particles come closer 
together and should be unity for large particle separation 
because a will not be influenced by far away particles. 


A useful analytic form is 


r: eo KV (4) 


er) Aoshi 


where V(r) is the Lennard-Jones potential and K a variational 
parameter. The next step is to express the expectation value 
of the ground state energy Ey for an assumed value of K in 
terms of a cluster expansion, truncate it after the first term 
and then minimize E> to get .a differential equation for the 
Single particle wave function do: This equation is solved 

and the ground state energy calculated as a function of K. 
Toelurns .utpthat 9, is approximated well by a spherically 
symmetric Gaussian (as is the case for the ground state of the 


harmonic oscillator) 
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where A is a variational parameter that minimizes Ey for a 
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given K. Finally, the ground state energy is minimized by 
varying K. The whole procedure is formally equivalent to a 
Hartree calculation with no correlations in which the Lennard- 


Jones potential is replaced by an effective potential 


2 
Vieg(t) = [V(r) + lini Piiatgen nie coos Se! ros 


This effective potential has a softer core than the Lennard- 
Jones potential. Using manybody techniques, Fredkin and 
Werthamer (1965) have shown how to calculate the phonon fre- 


quencies from these single particle exitations. They find for 


the eigenvalue equation at T = 0 K 
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The above equation is very similar to the eigenvalue equation 
2.9 of classical lattice dynamics with the exception that the 
interparticle potential has been replaced by an average with 


respect to the mean position probability of the two particles 


<0|V |O>. Pa at | 


eff 


From here on the equations of classical lattice dynamics 


(Leibfried and Ludwig 1961) apply. This is also a plausible 
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theoretical explanation why helium is in many respects a solid 
like any other which has long been established experimentally. 

The single particle picture has been refined in various 
ways but all numerical calculations were done for the bcc 
structure only. Hetherington et al (1967) extended the cluster 
expansion to include 3-body terms. They find that the ex- 
pansion converges. indeed sufficiently fast for a correlation 
of the form 2.13. Another expansion given by Brueckner and 
Frohberg (1965) however seems to converge less well. The 
ground state energy was also computed by Hansen and Levesque 
(1965) with a Monte Carlo Method using a different correlation 
function. Etters (1968) considers the dynamical behaviour of 
aj pair of particles in.the averaged field_of all. other. par- 
ticles. 

Another improvement especially in the volume dependence 
of the calculated quantities is expected from a better corre- 
lation function f. The usefulness of the particular choice 
2.13.for f.is mainly in.the rapid. convergence of the cluster 
expansion. Brueckner and Frohberg (1965) and recently 
Mullin et al. (1969) have tried different correlation functions 
obtained by numerically solving a differential equation for f. 
To find a better correlation function is difficult because the 


cluster expansion converges only for a restricted class of 
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functions f and not all variational parameters are suitable 
variational parameters. Massey and Woo (1968) use an effect- 
ive potential VoFe that includes liquid like correlations 
while in Meissners formulation (1968a, 1968b) both long and 
short range correlations are included. No numerical results 
have yet been published for sound velocities in hcp helium 
based on these improved correlation functions. 

The theories mentioned so far emphasize the single 
particle aspect of the problem and the phonons are constructed 
from the normal modes of response of the crystal to an ex- 
ternally applied disturbance. 

The socalled self consistent phonon theories (Hooton 
(1958), Boccara and, Sarma (1965), Ranninger (1965), Koehler 
(1967, 1968), Choquard (1967)) consider the phonons as the 
basic coordinates for describing the crystal. (thereby assuming 
that they exist in solid helium) and describe the motion of the 
individual atoms through the occupation number of the different 
phonon modes. While the two different approaches give similar 
results (Koehler (1967)) at T = 0 K, the collective theories 
seem to be more suitable for numerical calculations. In 
Koehlers work, this was achieved by expanding the crystal 
potential in Hermite polynomials (being the eigenfunctions 


of the harmonic oscillator) instead of a Taylor series. 
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The idea underlying the collective theories is the following 
(Koehler (1968)): 


The true crystal Hamiltonian A, 


is approximated by a harmonic Hamiltonian Hy 


where the entire set of spring constants G is to be determined 
variationally. Koehler finds that the choice of G that 


minimizes the energy is 
Geers Oe Pte 


and that the ground state eigenfunction of the harmonic 


Hamiltonian is 
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The similarity with the single particle picture is apparent 


and the main difference is that instead of averaging with 
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respect to single particle distributions, the average is here 
with respect to the phonon groundstate. The theory is self- 
consistent because the spring constants G depend on the 
eigenfunctions |0> which in turn depend on G. 

As in the single particle theory, a Jastrow factor 
has to be introduced to deal with correlations. 

Horner (1967) uses a diagramatic technique which in 
lowest order is equivalent to the selfconsistent phonon- 
approximation and also includes finite lifetime corrections, 
which affect the sound velocity via the Kramers-Kronig 
relations. This effect is also included in the papers by 
Ranninger (1965), Choquard (1967) and Werthamer (1967). 

Unfortunately, only two papers (Nosanow and Werthamer 
(1965), Gillis, Koehler and Werthamer (1968)) contain numerical 
values for sound velocities in hcp fore Nosanow and 
Werthamer's (1965) calculation is based on a single particle 
theory with a Jastrow factor of the form 2.13 and a Lennard- 
Jones inter-particle potential with © = 10.22 K and.o = 2.556 A. 
The cluster expansion of the groundstate energy is truncated 
after the two body terms and the energy estimated to be accurate 
within a few percent. 

Within the framework of the selfconsistent phonon 


theory Gillis, Koehler and Werthamer (1968) have calculated the 
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phononspectrum at two different densities and the sound 
velocities at the density of 19, 2cm>/mole; using the same data 
for the interparticle potential and the same Jastrow-factor 
as Nosanow and Werthamer. The internal energy is lowered by 
about 17 cal/mole in comparison with values given by the 
Single particle theory (Nosanow (1966) Hetherington et al. (1967)). 

Figure 2 shows the two theoretical results for the 
sound velocities in different crystal directions. They will 
be compared with experimental data in sections 4.1 to 4.3. 

Sound propagation in hexagonal crystals will be 
discussed in section 2.3 and the elastic constants obtained 
from the data of Gillis et al.will be used as an example to 
illustrate the concepts. 

In order to give an impression of the: performance of 
the other theories I show in figure 3 the internal energy E, 
at absolute zero for bcc He? for which most numerical calcu- 
lations have been done. The reason for this is that the cubic 
structure is simpler than the hexagonal, and only in He? does 
this structure exist at absolute zero and over a sizable 


volume range. The experimental points are from Pandorf and 


Edwards (1968). 
From E, (Vv) we obtain the pressure p 


p = -dE,/dv 2633 
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FIGURE 2 


Theoretical sound velocity vs. molar volume, according to 
Nosanow and Werthamer (1965) (solid lines), and Gillis et al, 
(1968) (circles, triangles and squares). y is the. angle 


between c-axis and direction of propagation. 
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the Bulk modulus B 


B= -V dp/dv = Vv. a*E /av" 2.24 


and an approximate expression for the longitudinal sound 


velocity Ye 


i) See Neaee (EL /Mhy: a*E/av* 2.25 


(V is the molar volume, M the atomic weight). 


2.3 Sound in hexagonal crystals 


Before measuring the anisotropy of the sound 
velocities in a crystal one would like to know something in 
principle about the velocity surface which consists of the 
three surfaces obtained by plotting in each direction the 
longitudinal and the two transverse velocities. They must 
have the same symmetry as the crystal and I will show that in 
addition they have to meet some other conditions. 

The surface is described completely by the | ae 
constants and the density . The elastic constants Ci5k1 
describe the stresses O55 in various directions in the crystal 


in terms of the strains ¢€ applied in various other directions: 
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FIGURE: 3 


Ground state energy E, vs. molar volume V in bec He? 


The curves. represent calculations by 


i. Deondon (g@54, p. 24) 

2. Nasanow (1966) 

3. Hetherington et al, (1967) 
4, Mullin et al. (1969) 

5. Horner (1967) 

6. Koehler (1967) 

7. Hansen and Levesque (1968) 


8. Iwamoto and Namaizawa (1966) 


The circles are measurements by Pandorf and Edwards (1968). 
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S44 is the force acting on plane i.in direction j.and the 


strain Ek is defined in terms of the displacement _u as 


du du 
sto Chk ate 
se ey (ite ne] fs 


ar k 


Equation 2.26 is nothing else than Hook's law. 
Knowing the forces on a small elementary cube allows us to 


write down the equation of motion for it 
2 
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el. =r eed = Ga. _=— C- 2228 
i ; Ax « 4k1 a Hg aa ere Gl 

This has a plane wave solution 
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which is nothing else than the leading term in the Fourier 
expansion 2.7. To get the sound velocities we then proceed 
in the same way as before by solving the secular equation 
for the desired wave vectors q. 


The number of independent elastic constants is greatly 
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reduced using symmetry and thermodynamic arguments and only 


5 remain for hexagonal crystals namely c c C and 


Tr 2” 


Cay (we have used the shorted notation of Voigt 


wR EA ape Fs 


iy cae gation Bx Be 


= Cc = C = Cay) and the 


sola bege: 12) “1033, “13! ©3333 337 - 9325 


matrix Cay for hexagonal symmetry is 


C17 S12 13 0 0 0 
Veo te his oO eRe? 
C13. °13 «$33 0 0 0 
ae 0 0 0 Cay 0 0 Bak 
0 0 0 0 Cay 0 
0 0 0 0 0 3 {11 7°12)) 


Several authors (Zener (1936), Gold (1950), Musgrave 
(1954)) have explicitly solved the secular equation for hexago- 
nal crystals and find that for this symmetry (and no other) 
the cubic equation factors into a quadratic and a linear 
equation for all directions and the sound velocity can there- 
fore be expressed explicitly as a function of the direction 
cosines of the wave vector g. Another result which is of 


decisive importance for the design of the experiment is that 


the velocity surface has a cylindrical symmetry around the 
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c-axis, which means that the sound velocities depend only 
on the angle y between wave vector and c-axis. 


The expressions for the sound velocities are: 


5 
Vv, = (c,/p) 
4 
=o Lae, 
Very Soee/e) 
1 
2 
vee (c, 7/0) 
with 
2 2 
c, = 4[(c +c ) sin y + (c +c ) cos y +¢ (y)] 
a me 4k 4 
c, =4(c - c ) sin?y +c cos* y e239 
C1 11 ray 
c, =4[(c +c ) sin*y+(c +c ) cos*y -¢(y)] 
Le ey 4 33 4 4 
Ma=(cy ec) sin yt? (ac -"c By ices? t yet 
Ll 4 33 rary 
aig 5 + 2 +7 : 
+ 2sin“y cos“y Bev Cre Co! (c CEs) (c cae ] 


For propagation along the c-axis (y = 0) 
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co = .¢c 

ds 1a 
c = $(c ee 
Ci 11 12 
Cc = Cc 

Ce 4 


A computer program for calculating velocities from elastic 
constants is found in the appendix. 

If the sound velocities are known as a function of 
the angle y, the elastic constants of the solid can be calcu- 
lated. The most convenient set of values are the velocities 
along the c-axis, perpendicular to the c-axis and for a 
direction in between (which is necessary to obtain Cy 3)- 

A section through a velocity surface based on theo- 
retical sound velocities published by Gillis et al. (1968) 
is shown in figure 4. The corresponding elastic constants 
are given in section 4.3. 

Some pecularities about sound propagation in crystals 
are worth mentioning here. One of them is that in general 
the direction of propagation of a sound beam does not coincide 
with the wave normal so that we have the situation shown in 


figure 5. The formula for the angle A; between beam direction 
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FIGURE 4 


Velocity surface of hcp He- according to Gillis et al. (1968) 
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and wave normal for the polarization i is given by Musgrave 
(954). 

Since it is rather lengthy I present it in the form 
of a computer program in the appendix. Figure 6 shows 
numerical values based on Gillis (1968) elastic constants. 

The deviations are surprisingly large and have to be 
considered in the design of a sample chamber. The chamber 
used by Vignos and Fairbank (1966) cuts off beams completely 
for which A» 10° thereby preventing the observation of long- 
itudinal sound for directions y in the range 20 ... 45°. This 
may partly explain that their data do not show the same large 
anisotropy as found in my measurements where deviations A up 
to 30° can be accomodated. Similar considerations apply to 
the chamber of Lipshultz and Lee (1965) which allows a devi- 
ation A£16° that cuts off almost all directions for the 
Soe branch. This property can be used to obtain elastic 
constants in bcc He? (Wanner (1970)). 

The three sound velocities in each direction have 
three different polarizations. In general, the polarizations 
are not "pure", that is the particle displacement for the 
longitudinal mode is not perpendicular to the wavefront, nor 


are the transverse motions in the plane of the wavefront. In 


hexagonal crystals however one transverse branch (usually 
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labelled ty) is always pure due to the rotational symmetry 

of the velocity surface. Its polarization is perpendicular 
to the c-axis and the propagation vector and it Heracles has 
the lowest anisotrophy of all branches. The other transverse 
branch t. and the longitudinal branch 1 deviate by an angle 

6 (figure 6) but the 3 polarization vectors are mutually 
perpendicular. 

The motion of a transducer surface will therefore be 
split into 3 components, and the transducer will in general 
excite all three branches with an amplitude proportional to 
the corresponding component along the polarization vector. 

It is difficult to observe the transverse components 
when exciting with a longitudinal transducer because the first 
transverse echos often coincide with the second longitudinal, 
but longitudinal sound coming from transverse transducers 


has occasionally been seen both by us and Lee et al (1969). 


2.4 Related Quantities 


If the transverse and longitudinal velocities are 
known for a sufficient number of directions, the elastic 
constants can be determined. From known elastic constants 


several other quantities such as the Debye temperature 0, 
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FIGURE 5 


Propagation of longitudinal and transverse sound beams in 


hexagonal crystals (schematic). 
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FIGURE 6 


Deviation A; of beam direction and deviation 6 of polarization 
vector from wave normal vs. orientation y. The curves are 
calculated from elastic constants obtained from the work of 


Giifis et al.'(1968).. 
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at 0 K and the compressibility can be calculated and compared 


with thermal measurements. 


Debye temperature Oe 
The Debye temperature at 0 K is (Alers (1965)) 
i 
eg ea kaa 
O k \4nv ‘ 


where the Debye velocity v is 


C -3 =3. -3 =1/3 
V = (SF (vy SUNS Day) aQ/12T) 
h = Planck constant 
k = Boltzmann constant 
N = Avogadro constant 
V = molar volume 


2<34 


eget) 


The integration over the solid angle 2 is particularly 


easy for hexagonal crystals, where the only integration 


variable is the angle y between wave vector and c-axis. 


A 


computer program for 0, using the elastic constants as input 


is given in the appendix. The Debye temperature so obtained 


should be identical with the calorimetric Debye temperature 


measured at temperatures below 0/90. 
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Compressibility, isothermal and adiabatic elastic constants 


For hexagonal crystals the isothermal compressibility 
K = -(dV/dp)/V in terms of the isothermal elastic constants is 


(Nye (1946), p. 146) 


Cc +c - 4c + 2c 
om Li is, 13 35 : 2 a6 
Spates, 'qcxoGase vere 


To compare the compressibility obtained by thermo- 
dynamic measurements with the compressibility calculated from 
sound velocity data one has to keep in mind that the former 
are static measurements giving the isothermal compressibility 
while the sound velocities give adiabatic elastic constants 

Ss 


and compressibility. The difference Cras cee depends on the 


particular elastic constant and is of the order (Bhatia (1967) 


ike Sok eee A ale 2037 


(C_, C.. heat capacities at constant pressure, volume). 


Further comments on this point are found in section 4.4. 
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Phonon spectrum 


The sound velocity is equal to the phase velocity of 


phonons 


Max 


(w = phonon energy, q = wave vector, Gage = MAYS VEC lor ar 

zone boundary) in the limit of long wavelengths Leibfried 

and Ludwig (1961, p. 426) argue that the above definition 

may give isothermal sound velocities rather than adiabatic 
velocities, but leave the question open. In neutron. scat- 
tering experiments, the excited phonons have usually rather 
high frequencies compared with ultrasonic frequencies and 
elastic constants calculated from 2.38 are therefore isothermal 
constants (see section 4.4). I have purposely avoided the 


notation "q. > 0", since this, taken literally, would give 


adiabatic constants. 


Second sound 


Assuming an isotropic material, the velocity of second 


sound is (Ackerman and Guyer (1968) ) 
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Cauchy relations 


tt ’every atom in’ the crystal “is a°centre of inversion 
Symmetry and if there are only central forces between a pair 
of atoms the elastic constants obey relations known as Cauchy 
relations (Born and Huang (1962)). For hexagonal crystals 


these are 


Q 
I 
WwW 
Q 


In a hexagonal close packed crystal the condition of inversion 


symmetry is not satisfied, but Leibfried (1955) shows that 


this affects only the relation C11 = 3C4 5: For the case of 
next nearest neighbour central forces C13 = Cyy remains 
valid. 


2.5 Elastic constants from unoriented single crystals 


in the cubic phase 


Finding the orientation of the crystals with the 
birefringence technique is possible only for the hexagonal 
phase. The orientation of cubic crystals 1s done either with 


neutron diffraction or with x-rays. It is believed that 
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neutron diffraction is not a suitable method for hee because 
of the large capture cross section for this isotope and x-ray 
scattering has so far been applied only to samples of a few 
millimeter diameter, because the small x-ray scattering cross 
section of helium requires thin container walls. On the 
other: hand bcc fer is of great experimental interest. The 
large deviation of the specific heat from the Debye theory 
at low temperatures is as yet unexplained and sound velocity 
measurements at low temperatures would at least establish 
the phonon contribution to the specific heat. Froma 
theoretical point of view, this phase and isotope have been 
studied most extensively. I have therefore investigated the 
possibility of finding the elastic constants of a cubic 
crystal using ultrasonic means alone. Assuming that the 
crystals do not grow in preferred directions, that many 
crystals can be grown and a substantial fraction of them are 
single crystals, this should indeed be possible. 

The method consists of measuring the longitudinal 


velocity Vi and at least one transverse velocity v, along the 


ie 
same path and then plotting them on a (Vir Fela graph. For 
cubic single crystals the points are restricted to the two 


Shaded areas indicated in‘’figure 7. . These areas are 


topologically identical for all cubic crystals and the co- 
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ordinates of their prominent features, namely of the corners 
and extrema can be expressed in terms of the elastic constants. 
For most cubic elements the figure lies in the (vir fa plane 
as drawn, but for some binary compounds it may be rotated by 
180°. The calculations (Nosanow and Werthamer (1965), de 
Wette et ale(1967)) show that in this respect helium is 
Similar to most of the elemental crystals, but should the 
figure turn out to be upside down the extension of the analysis 
would be straight forward. 

Sound velocities in polycrystals will increase the 
areas and fill in the forbidden zone in the center but the 
extremas and corners will not be obscured. With a sufficient 


number of (Vi, v,)- points it should be possible to determine 


ms 
their coordinates and calculate the 3 elastic constants from 
them. Five linear combinations of Cys C15 and Cay are 
obtained which allows a check for consistency. This method 
works also for hexagonal materials, even better than for cubic 


Since the (v Ved points are restricted to ly on two curves 


VW 
rather than in two areas. 

Once the elastic constants are obtained it is in 
principle also possible to find the orientation of. the 


crystal from sound velocity measurements alone (Green and 


Henneke (1967)). 
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FIGURE.7 


Transverse vs. longitudinal sound velocity in bcc He wl based 


on elastic constants from de Wette et al. (1967). 
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3. EXPERIMENT 


3.1 Low temperature apparatus 
SpLan Cryostat 


Helium solidifies only under pressure. The high 
pressure cell in which the helium crystals are grown must 
therefore contain transducers to measure the sound velocity 
and a window through which the orientation of the crystal can 
be determined. In order to grow and anneal a crystal, the 
temperature of the cell must be controlled during a long time 
within close tolerances. To meet these design parameters, 

a Special cryostat was built. 

Crystals are grown by keeping the temperature at the 
top of the pressure cell slightly above the melting temperature 
with an electric heater and by slowly decreasing the bottom 
temperature by pumping on a helium pot on which the cell is 
mounted. The cell is in vacuum and the high pressure fill 
line is vacuum jacketed to prevent it from blocking. Figure 
8 shows the layout of the pumping system. An important part 
of it is a pressure regulator (Walker (1959)) that can be 


connected either to the pumpline for the dewar or to the 
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FIGURE & 


Pumping system 


1. Wallace - Tiernan precision pressure gage 
2. Walker pressure regulator 

3. Needle valve 

4. Ballast volume 

5. High vacuum gage 

6. - Drtfusion pump 


7. Exchange gas container 


¥. Vent valve 


R Connection to helium gas recovery system 
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pumpline for the pot. Usually the helium in the dewar was 

at atmospheric pressure and 4 K, the regulator. connected to 
the pot pump line and the dewar pump used in parallel with 
the pot pump. The regulator holds temperatures within a few 
mK in the temperature range 1.4 ... 4 K. It can very easily 
be adapted to automatic crystal growth (Ackerman (1967)) by 
including a low flow calibrated needle valve through which 
the reference pressure in the ballast volume is slowly pumped 
down and the pot temperature therefore slowly decreased. 

One stilling of .the. 230 om? pot lasts about 16 hours 

at 2 K and it can be refilled through a needle valve in about 
10 min. while pumping.on it. In this way it is.possible.to 
keep the temperature below 2 K for an arbitrary long time. 
A rough estimate of the pot temperature is obtained from the 
vapour pressure which is read on a bourdon gage (Model FA 
145, Wallace & Tiernan, Belleville, N.J.) in the pump line. 
The lowest temperature obtained was 1.23 K. The pot is sus- 
pended on two stainless steel tubes, one of which is. the 
pumpline and the other serves as a vapour pressure bulb for 
calibration of a germanium thermometer which is attached to 
the pot near that bulb (Figure 9). 

Optical access into the copper can is through a 


vertical 13 mm tube with a room temperature window at the top 
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of the cryostat. This is a very easy solution from the point 
of view of cryostat construction but requires a window or 
prism inside the cryostat and one on top of it to get a 
horizontal light beam. This complicates the optical problems. 
A nickel vane inside the vacuum container and at low temper- 
ature can be moved from the outside with a magnet to block 
the optical path so that no room temperature radiation enters. 
However neither the temperature control of the pot nor the 
growth of helium crystals was in any way affected by opening 
or closing this vane. 

Shortly after completing the measurements on longi- 
tudinal sound velocity I had the unfortunate experience of 
observing pressure oscillation in the helium gas above the 
liquid level of the dewar similar to the oscillation in the 
familiar dipstick level indicator. They increased the 
evaporation rate from the previous 140 cm?/hr to about 2500 
om?/hr and prevented further experiments. I have tried to 
suppress these oscillations by a styrofoam plug about half 
way between the liquid helium surface and the top of the 
dewar, by filling the brass head of the cryostat with styro- 
foam or. aluminum foil, by connecting a 10 litre ballast volume 
to the dewar, by disconnecting the dewar from the recovery 


line and venting it to atmospheric pressure, and by using a 
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FIGURE 9 


Cryostat 

1. Copper can, enclosing high vacuum 

2. Vacuum jacketed high pressure capillary 
3. High pressure sample chamber 

4. Connector to sound transducer 

5. Germanium thermometer 

6. Pot, filled with pumped liquid helium 
4. Laser beam 

8. Allen-Bradley thermometer 

9. Manganin heater 

£07 Mirror 

11. Adjustable mirror support 
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slightly longer’ dewar (78 cm instead of 73 cm), but with little 
success. The oscillation would sometime be absent only to 
appear again spontaneously after the next filling of the dewar. 
The inside of the dewar is 9 cm and it is connected to a brass 
head of 12 cm inSide diameter and 13 cm length through. a 
Corning flange (Pyrex conical pipe system, Corning Glass works, 
Corning N. Y.). The round copper can inside the dewar measures 


Pad Xeb267em.. 


3.12 The high pressure cell 


The aim of the experiment was initially to measure 
not the velocity but the absorption of sound and I spent there- 
fore much time in developing a suitable support for the guartz 
transducers. With this support I have observed up to 72 sound 
echos in liquid helium and up to 200 in solid helium. In the 
first bomb that had an optical window (figure 10) I retained 
this support and did some preliminary measurements at a pressure 
of 58 and 128 bar with it. This chamber hangs freely from 
the pressure capillary inside a copper can which is submerged 
in a pumped bath of liquid helium. Crystals were grown by 
cooling the chamber with exchange gas and keeping the capillary 
from blocking with a manganin heater wire that runs inside it. 


Not surprisingly, the success rate for making single crystals 
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FIGURE 10 


Sample chamber 1 


1. High pressure capillary 


2. Lens ring seal 
oe Sprang 
4. Sound transducer assembly 


5. Quartz transducer 
oe.  2ndium O-ring 


7. Quartz window 
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was quite low and led me to design an improved version which 
is described below. The crystals in the preliminary chamber 
were observed from the top of the cryostat through a window 
in the bottom of the cell using a system of two mirrors 
attached to the cell. 

In the improved sample chamber (figure 11), crystals 
were grown with a vertical temperature gradient and I made 
the bottom of the cell from copper which has a high- heat 
conductivity and can carry away the heat of solidification 
to the helium pot. The cell itself is made from a 2% beryllium 
copper alloy which has high tensile strength, low heat con- 
ductivity (compared with copper) which makes it easy to pro- 
duce a temperature gradient, and is easy to machine, polish 
and solder. After machining it is hardened for 1 1/2 hrs at 
B15 ic. 

The special shape of the bottom should induce nucle- 
ation at the coldest spot A. One of the nuclei will then out- 
grow the others as the solid liquid interfaces moves. up and 
act as a single crystal seed for the larger crystal higher up. 
I have however no proof that this arrangement really increases 
the probability of single crystal growth. It is around 50% 
with the growth technique described later. 


No obstacles should be in the path of the growing 
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FIGURE ll 


Sample chamber 2 


1. High pressure capillary 

fae tendium: O-ring 

3. Electrical connection to backside of sound transducer 
4. Stycast epoxy for electrical insulation 

be © Solid: Heirum 

Oe OOLt, solder 

7. Helium pot 

8. Sound transducer 


9. Window 
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crystal and Dr. J.E. Vos (Delft) has told us that he sometimes 
observed pinning of grain boundaries to discontinuities in 
the wall of the container. My cell is a tube with window and 
transducer mounted in holes drilled from the side into it. 
The sound transducers are thin quartz plates of 12.5 mm 
diameter plated on both sides with a layer of aluminum as 
electrodes. The transducer serves also as a mirror to reflect 
the light beam that is used for orientation determination out 
of the cell again. Therefore the cell requires only one window 
which facilitates corrections for stress birefringence, and 
this window is at the same time a reflector for the sound beam. 
The advantage is that light and sound beam coincide in 
direction and position and sample the same part of the crystal. 
I used the greased indium O-ring technique described 
by Vos and Kingma (1967) to seal the fused silica window to 
the bomb. The O-rings had initially a square cross section 
of 0.7 mm and an inside diameter of 15 mm. After compression 
the thickness was approximately 0.1mm. The sealing surface 
on the bomb was hand lapped and polished to a mirror finish. 
This seal withstands repeated temperature and pressure cycling 
but fails eventually. The same technique is used to seal the 
backing plates for the transducers. to the cell. The center 


piece of the plate is electrically insulated from the cell 
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and in contact with one electrode of the transducer. The 
other electrode facing the helium is. grounded to the bomb. 
This centerpiece is epoxied into place with Stycast 2850 GT 
(Emerson & Cumming Inc., Canton, Mass.). The tapering should 
prevent back reflection of sound. The surface against which 
the transducer is pressed with to washer springs is lapped 

and polished until a check for flatness with an interferometer 
plate in sodium light shows less than one fringe over the 
whole surface. This is more than adequate for a sound wave 
length of approximately 0.1 mm. 

The parallelism between window and transducer is not 
very critical for velocity measurements. Coarse adjustment 
was achieved by shining a laser beam onto the transducer and 
compressing the indium O-ring until the reflection from the 
transducer coincides with the reflection from the inside face 
of the window. With careful tightening of the screws inter- 
ference fringes are observed and circularized. I estimate 


window and transducer to be parallel within about Tae radians. 


3.13 Temperature measurement and control 


I determined the temperature of the helium pot by 
measuring the resistance of a germanium thermometer (Solitron 


Devices, Riviera Beach, Fla.) with an isolating potential 
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FLGURE*12 


Electronic temperature regulator 


Rl. Allen-Bradley thermometer, 3.5 k at T= 2.5 K 

R2 Manganin heater 

R3 Dummy heater 

Sl Thermometer current on/off 

S2 Time constant selector 

S3 Heater current on/off 

S4 left side: measure Rl, manual heater current regulation 
right side: automatic heater current regulation 


105 Operational amplifier, model 105 (Analog Devices) 
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comparator circuit (Dauphinee and Mooser (1955) Dauphinee and 
Woods (1955), Dauphinee and Preston-Thomas (1958)). The resistance 
at 2 K is about 1000 Ohm, the measuring current to Amp. and 
the precision about 0.1 mK. The thermometer is calibrated 
against the vapour pressure of He? using the Te 9 scale 
(Brickwedde et al. (1960)). For unknown reasons, but most 

likely because I did not wait long enough for thermal equi- 
librium, the calibration points show a scatter of + 5 mK from 


a smooth curve. For interpolation a polynom of degree 6 is 


fitted by computer to the points using the method of least 


squares: 
: 1 
Logan. = slew As 10g. RB) cide 
oir 
ie SO) OO a 41.815, A, = Oe Lats A, = Se2OL85 
a ee 0.42763, A, = 20. L94104 ait 0.01638 119 


The temperature at the top of the pressure cell is 
found from the resistance of an Allen-Bradley radio resistor 
with a room temperature resistance of 100 Ohm. It is cali- 
brated against the germanium thermometer using exchange gas 
to establish thermal equilibrium between the two thermometers. 


This thermometer serves also as the sensing element in an 
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electronic circuit that keeps the top temperature constant 
and above the melting temperature to keep the pressure capillary 
from blocking and provides the necessary temperature gradient 
for crystal growth. It forms one arm of a Wheatstone bridge. 
The out of balance dc-signal is amplified and integrated by 
an operational amplifier with a time constant that can be 
selected between 0.01 and 250 sec. and then amplified further 
by a power amplifier which controls the current through a 
100 Ohm manganin heater. This feedback loop regulates the 
top temperature around 2 K to within a few mK when cooling 
the bottom of the cell to 1.3 K and has a drift under stable 
conditions of less than 1 mK over a period of 8 hrs. Used 
together with a dummy heater and a meter for the heater 
current it also serves as a thermometer. I would like to 
thank Dr. J.S. Rogers for help in the design of it. 

It turned out that too much heat is flowing down the 
Capillary (stainless stell, OD = 2.5 mm) from room temperature. 
A copper wire clamped to it and anchored at the pot is enough 


however to reduce the top temperature. 


3.14 High pressure gas handling system. 


Mezhov-Deglin (1965), who was the first to successfully 


grow single crystals of helium, stresses the importance of 
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purity of the initial gas and reports that best results are 
obtained by cleaning it with a leak that allows passage of 
superfluid helium only. On the other hand Ackerman and Guyer 
(1968) have observed second sound in crystals prepared from 
gas obtained from evaporating liquid helium. The sound 
velocity is probably much less affected by impurities than 
the sound absorption and it has to be only sufficiently high 
to allow the growth of single crystals. 

The gas handling and cleaning system is shown in 
figure 13. The sample chamber is evacuated and flushed several 
times at room temperature with cleaned gas and then cooled 
with exchange gas under constant pressure of about 10 bar to 
77K. Then it is slowly pressurized to either 36 or 58 bar 
and cooled further to 4 K, thereby sucking gas through the 
U-tube that is dipped into a storage dewar. Cooling from 77 
K to 4 K takes approximately one hour which should be suffi- 
ciently slow to liquify all helium that passes through the 
trap. A storage cylinder is used as ballast volume to maintain 
constant pressure during the growth as more helium passes into 
the cell. Pressure stability is therefore mainly determined 
by the amount of room temperature fluctuation and always less 
than 0.1 bar. I read the pressure on a Bourdon gage (Heise 
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S215 Crystal Growth 


Crystals were grown at constant pressure and with a 
constant growth rate of 2 cm/hr using a technique described 
by Ackerman and Guyer (1968). I would like to thank Dr. 
Ackerman for telling me about it prior to publication. The 
main idea is to keep the fluid helium slightly above the 
melting temperature (about 20 mK) and to lower the temperature 
of the solid helium in such a way that the solid-liquid inter- 
face progresses upwards at constant velocity. Since the heat 
conductivity of the solid is several orders of magnitude 
higher than the liquid conductivity it can be assumed that all 
the heat of melting is carried away through the solid (the 
walls have a similar or smaller conductivity as solid He and 
are ignored). From the known temperature dependent conductiv- 
ity (Bertman et al.1966) the temperature at the bottom can be 
calculated as a function of time for a given growth wares 

Liquid helium I has a temperature conductivity 
k/o.C (k = heat conductivity, C.= heat capacity per. gram) of 
about ie" om? gece so that it takes about ido sec to reach 
thermal equilibrium in a volume of l com>. Convective flow is 


absent because of the very small thermal expansion of liquid 
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FIGURE 13 


High pressure gas handling system 
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helium. Usually I made a new crystal after melting the previous 
and waiting for one or two hours for equilibrium. I have tried 
to observe the influence of that waiting time on crystal quality 
but find no conclusive results, probably because it depends 
also on many other factors. However after completing the 
experiments I calculated the average waiting time for the 
runs that yielded single crystals and find 90 min while the 
time for polycrystals was 65 min. Ackerman (1968) waits 6 to 
10 hrs with a chamber of similar size. 

I keep the sound switched on before and during the 
growth to monitor the progress of the crystal. Since at 36 
bar the sound velocity is temperature dependent near the 
melting line due to the vicinity of-the A -transition, this 
also allows to check for temperature equilibrium in the liquid. 
After a growth time of 60 min the liquid echo decreases in 
amplitude in a fluctuating way but the time delay remains the 
same. This shows that the pressure inside the cell has remained 
constant. After 100 mina small, fast echo from the solid is 
visible and grows at the expense of the liquid echo or, quite 
often, the liquid echo disappears completely before a new 
echo can be seen. After about 150 min the solid echo is stable 
in amplitude, or grows very slowly (figure 20). 

Another indication shows that solid helium is in 


contact with the transducer: The transducer oscillates long 
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FIGURE 14 


Temperature at bottom of sample chamber vs. 


constant.growth rate of 2. cm/hr. 
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after the 2 usec electrical high frequency pulse is stopped 
because of the high Q-factor. The decay time depends on the 
acoustic impedance of the load and is about 53. usec for liguid 
and 35... 39 usec for solid helium, just about inversely 
proportional to the acoustic impedance. It should in principle 
be possible to measure the sound velocity at the surface via 
this effect and then compare it with the bulk sound velocity 
to. determine crystal quality. 

After the whole bomb is filled with solid, I. reduce 
the heater power to zero and simultaneously raise the bottom 
temperature to within 0.2 K of the melting temperature to get 
equilibrium. Between 2 and 4 hrs after beginning of growth 
the velocity usually changes, mostly to higher values and 
seldom more than 1% in good crystals and then reaches a steady 
value. I assume that during this time recristallization takes 
place and the growth of large single crystals which eventually 
fill the whole cell. Annealing for a longer time or closer to 
the melting temperature does not increase the echos any further 
and often decreases them. 

I have occasionally attempted to measure the temperature 
dependence of the sound absorption. In figure 15 I plot the 
logarithm of the intensity of the first echo versus temperature. 


The experiment extended over 3 days. While this measurement 
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FIGURE 15 


Intensity of longitudinal sound echo under repeated 
temperature cycling. The crystal was grown at a pressure 


of 36.2 bar, the sound frequency was 5 MHz, the total path 


length 4.57 cm. 
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is wildly irreproducible one can see that repeated temperature 
cycling from 1.9 K to 1.3 K and up again improves. the crystal, 
more so than a stay at- constant temperature. The cause may 

be that the temperature gradients which are invariably pro- 
duced when cooling and warming offer a preferred direction 

for the movement of the grain boundaries, and it appears 
therefore that annealing under a gradient yields better 
crystals. I would like to thank Len Vienneau for an illu- 
mMinating discussion on these points. 

Lee, Heybey and Crepeau have recently undertaken the 
same measurements, at a lower density, with crystals grown 
from superfluid helium by pressurization. They report (Heybey 
and Lee (1967)) that single crystals can be grown in less than 
five minutes. . The reason is that the heat of solidification 
is carried away much more rapidly because of the high heat 
conductivity of the superfluid. Since many crystals are 
required for a determination of the orientation dependence of 
the sound velocity, the importance of quick crystal growth 
can hardly be over estimated. The disadvantage of this method 


is that only a very limited density range can be studied. 
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See LODticCS 


3.2.1 Birefringence 


4 was determined with 


The crystal structure of solid He 
x-rays (Keesom and Taconis, (1938), Mills and Schuch, (1961), 
Mills and Schuch, (1962), Schuch and Mills, (1962), Schuch. and 
Mills (1963) and with neutron scattering (Henshaw (1958), 
Lipschultz et al. (1967), Minkiewicz et.al.(1968) Brun et al. 
(1968), to be hexagonal close packed over the greater part of 
the phase diagram. 

Because of the hexagonal symmetry helium is uniaxially 
birefringent which means that the index of refraction is 
different for light polarized parallel (= n2) and perpendicular 

= no) to the optic axis. The optic axis of a plate is the 
projection of the c-axis on the plate (figure 16). Plane 
polarized light will be split into two components and after 


passage through the plate they will have a phase difference $6 


d is the thickness of the plate, A the light wavelength. The 
resulting light.is in- general elliptically polarized. I -will 
in. the rémainder. of-this- chapter call any optical element that 
introduces a phase difference between two perpendicular 


components a phase plate and » the retardation of the plate. 
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The difference in the refractive indices, ne - oe 
depends on the direction of the lightbeam with respect to the 
crystal axis and in hexagonal crystals is a function of the 


angle y between beam and c-axis only 


Zz a R ee 
ne ae ae (n, n,) yo eae - ee | 


Strictly speaking this relation holds only for 
1) ale << 1 which is certainly true for helium (n-n, crs Me 
By measuring the phase difference one can therefore determine 


the position of the c-axis in space and the angle y if the 


birefringence is known: 


oh 
‘ef ad) 2 3.4 
a GE 


ee Nn, was measured by Vos et al.(1967 b,c) for densities 


between 11 and 21 cm?/mole and by Heybey and Lee (1967) at 


ak one Pore (see Table l.). 


TABLE l 
p n, 7 7, Ref. 
2 >eatm 2. GortabOndt 2% 1076 Heybey and Lee (1967) 
26 2e8cthr lad Vos et al. (1967 b) 
30 2Ze764th 90408 Vos et al.(1967 c) 
72 2 due &, (One O:7 Vos. et al- (1967 c) 
140 Aes OO ys Vos et al. (1967 c) 
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I want to emphasize at this point that it is a 
coincidence that we can use birefringence to orientate the 
crystals for ultrasonic work, since birefringence does not 
give us the position of the a-axes. Fortunately this is not 
required because of the cylindrical symmetry of the velocity 
surface in hexagonal materials. Cubic materials are neither 
birefringent nor do they have a simple velocity surface. 

One can measure the retardation in principle with one 
of the many methods suggested in textbooks of optics. These 
methods become however inaccurate or unreliable in the 
presence of straylight of which we have plenty due to the two 
windows between the helium crystal and the analyzing equipment. 
I use therefore a similar procedure as Vos et al (1967 b). A 
plane polarized beam of arbitrary polarization vector is 
sent through the crystal. The resulting light is elliptically 
polarized. With a quarter wave plate and an analyzer I then 
determine the two parameters that characterize it, namely 
the direction a of the major axis of the ellipse (measured 
from some arbitrary zero direction fixed to the laboratory 
frame) and the ratio tan p of the minor and major axis. I 
plot these parameters for about 25 different directions of the 
initial polarization vector, each approximately Kaye apart. 


The way a and 0 change with initial polarization obviously 
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Ellipticity and axis direction of elliptically polarized 
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depend on the properties of the plate, that is its retardation 
@¢ and the direction G of its optical axis and they can be 
inferred from the p(a) - plot (figure 7). 

I now return to a more detailed description of the 
individual steps in this procedure. The first is the 
measurement of p and a. To do this I place a quarter wave 
plate, which is just a phase plate with 6 = 90° with its axis 
parallel to the axis of the ellipse and remove thereby the 
90° phase shift between the two perpendicular components. 
After the quarter wave plate the elliptically polarized light 
is plane polarized. With an analyzer I then determine the 
Girection Of polarization by setting it. perpendicular ,to_1t 
for complete extinction. Of course I do not know initially 
Woere the-axis of the ellipse is.and have to. find by trial 
and error the position of quarter wave plate (=m, ) and analyzer 
(=n), that results in complete extinction of the light. I 


then get 


This assumes that both quarterwave plate and analyzer are 


properly mounted, that is the zero reading on their scale 
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coincides exactly with their optical axis, but usually this 
is not so. Fortunately I get complete extinction for two 
different but symmetric settings of plate (m5, m. ) and 
analyzer (ny, no) and can use them to calculate and eliminate 


that error 
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The next step is the calculation of the retardation 
@ and the axis direction G of the phase plate from the o(a)- 
plot. For this I make use of formulas given in Born (1965, 
p. 27) that relate the shape and position of. an ellipse to the 
phase difference and amplitude of two perpendicular oscillations 


and find for po(a) 


ae 2 
$1729 fo, tec} amsin —- 2 (a-G) Ss7 
cos $¢+tan’ 2(a-G) 


One can see that p = 0 for a = G and G can therefore be read 


@irectiy from the plot.” To-find the retardation ¢ we note 


that 
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Lop g 


A better way is to use the computer and find the 
parameters ¢ and G in equation 3.7 that give the best fit 
to the measured curve p(a). The problem has no analytic 
solution like many non linear least square fits and requires 
an iterative procedure. 

The optical observation in polarized light allows also 
to decide whether I have grown a single crystal or not. Un- 
fortunately I am too far away from the crystal, and the various 
windows and mirrors distort the image so much that I could 
not see the crystal directly. However, if more than one 
crystal is oPee che it is impossible to find a setting of 
quarter wave plate and analyzer that results in complete ex- 
tinction of the beam. Extinction could only be accomplished 
for one crystal and this setting would not block light of 
another polarization. The sensitivity of the criterion is 
given by the amount of stray light which I estimate at a few 
percent, so that stray crystals which cover less than a few 
percent of the field of view will go undetected. 


Another indication of single crystals is the fact 
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that the measured o(a) -curves obey the relationship 3.7 to 
within the accuracy of measurement of p and oa. A large 
scatter is usually due to difficulties in locating the 
extinction setting and means therefore that more than one 
crystal is present. For single crystals the rms deviations 
from the curve are of the same order (2... ) as in the 


measurements with liquid helium in place of the solid. 


3.2.2 Optical apparatus 


Figure 18 is a diagram of the optical apparatus used 
to measure the retardation. The light of the Spectra Physics 
Laser (Mod. 130) is monochromatic with a wavelength of 
632.8 nm, well collimated (beam spread 7 Yom rad) and plane 
polarized. The output power is 1 mW. Since the direction of 
polarization is fixed, a polarization rotator was built 
consisting of a quarter wave plate and a polarizer. The plate 
has its axis 45° from the polarization direction of the laser 
and produces therefore circular polarized light which is then 

plane polarized by a plastic sheet polarizer. The beam is 
then reflected into the cryostat by a mirror that can be 


rotated around a vertical axis and tilted around a horizontal 


axis with adjusting screws. The window is a glass disc of 
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FIGURE 18 


Optics for the birefringence measurement (schematic) 
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12 mm aperture pressed against a rubber O-ring by the atmos- 
pheric pressure. Inside the cryostat is another mirror that 
can be turned around the same axis as the mirror on top and 
also adjusted around a horizontal axis. The beam then enters 
the pressure cell through a 12 mm thick window of fused 
Silica, traverses the cell and gets reflected on the alummized 
surface of the sound transducer. The incoming and outgoing 
beam are at a slight angle (5x10 *rad) so that the outgoing 
beam just misses the laser and can be analyzed with a quarter 
wave plate and an analyzer. The extinction of the light beam 
is observed on a ground glass. Since the beam passes the 
crystal twice, its effective length in the crystal is doubled 
to 4.57 cm and so is the total phaseshift $9. 

As mentioned before the simple mechanical arrangement 
leads to a number of optical problems. The reflection at the 
surface of a metallic mirror introduces a phase shift between 
the components of the electric vector that are parallel and 
perpendicular: to the plane of incidence and the reflection 
coefficient for the two components are not exactly the same. 
I have tried to use totally reflecting prisms instead, because 
although they give also a phase difference for the two compo- 
nents at least the reflection coefficient is the same for 


both. I found the prisms unsuitable however because it was 
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difficult to keep the lower prism at constant and uniform 
temperature and it showed stress birefringence varying with 
time and temperature due to nonuniform thermal expansion of 
the glass. 

If we are only interested in the state of polarization 
of a light beam and not its direction we can represent a 
mMirror*through a phase plate. Its axis is in. the plane of 
incidence and its retardation is equal to the phase difference 
@ for reflection of parallel and perpendicular component and 
can be calculated from electromagnetic theory. The reason is 
that it does not matter for a mathematical description of the 
beam by which mechanism the shift has been produced in an 
element. (I have assumed here that the difference in re- 
flection coefficient between the two components can be neg- 
lected. It is 7% for aluminum mirrors and its influence will 
be discussed later). 

The problem is now how to measure the birefringence of 
a plate that can only be observed through a series of phase 
poaces, “two Of° it beingsthe? ’mirror’-platessand one the. high 
pressure window which might show stress induced birefringence. 

Vos et al. (1967 b) have solved that problem when only 
one plate is present, namely the birefringent window. The 


solution is based on a matrix calculus invented by Jones 
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(1941 a,b,c) in which optical elements are represented by 
matrices and a series of elements by matrix products. It 
involves measuring the retardation ~ and the axis direction 
p of the combination of helium crystal and window and the 
retardation 8 and axis direction b of the empty cell. The 
retardation ¢ and axis direction g of the helium crystal 


alone is then 


‘riage Reypreinten gen jalle se vaes cine us 
cos 5¢ = cos2(b p).sin5)-sin5B + coss.cos56 


Sin2(b-g) = sing) : sin2(b-p)/sing 0) 


This solution is easily adapted to the case of 
compound windows. Jones (1941 b) shows that a series of phase 
plates can be represented mathematically by one single phase 
plate and a rotation of the coordinate system. The rotation 
affects only the axis direction g of the crystal, in which 
we are not interested anyway, but otherwise the same formula 
as for the single window case apply, only £8 and b are now 
the retardation and axis direction of the whole optical 
system but without the helium crystal, or with liquid helium 


mich is not birefringent. 
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3.2.3 Errors in the orientation measurement 
pees ees oS Bee ES OES NCation.measurement 


The assumption made above that for metallic reflection 
the reflection coefficient is identical for components parallel 
and perpendicular to the incident plane needs some justifi- 
cation. It is very difficult to treat the whole optical 
system mathematically without making that assumption and I 
will therefore only examine what happens to the o(a) -curve 
of a single metallic mirror and an angle of incidence of iba 
Then I will verify experimentally that eqation 3.9 can indeed 
be used for a train of mirrors. For an aluminum mirror at 45° 
incidence the retardation is 168° and the ratio of the re- 
flection coefficients is 0.927. 

The axis direction of the equivalent phase plate 
remains unchanged, since this involves only the reflection of 
one component, and the reflection coefficient for the other 
does not matter. If both incident components are equal, that 
is fora =Gt 45° and one of them is reduced by 7% more than 
ines other through. the reflection,: this is.equivalent to.a 


1 0.927) 


rotation of the coordinate system by an amount (fan) 
= 45 in other words a distortion of the a -axis by a maximum 
of 2°. This is within the accuracy of measuring a. 


For the experimental verification I put a mica phase 


plate at room temperature in place of the helium crystal and 
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measure its phase difference through the cryostat. Giving it 
3 different orientations gave the same results within 1° and 
agreed with a direct measurement to within 1°. The retarda- 
tion: of- that ‘plate ‘was 28°. With a 67° -plate, consistency 
was again 1° and the deviation from the direct measurement 2°. 
The total retardation B of the empty cryostat was 32° at 

room temperature. 

During the actual run 8 changed to 80° and the empty 
cryostat had to be measured at frequent intervals. This could 
be due to formation of an oxide layer on the mirror or a 
change of the complex refractive index of the bottom mirror 
with temperature, and an uncertainty of 5° for @ is probably 
more realistic. The corresponding error in the orientation y 


O 


is greatest for 8 = 0 and 6 = 90° since dy/d¢ is infinite 


there. This is an unfortunate aspect of the birefringence 
method since these are just the direction of prime interest. 


Around y = 10° ana Y= 80° the uncertainty is of the order of 


no” | going down to + 5° at y= 45°, 


The error due to the uncertainty of crystal thickness 


is negligible compared with this, but the uncertainty in the 


birefringence LS oe dearest of about 3% contributes 1° at Y= 45°, 
O 
going up to 7° at yere= SY OP 
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Sao. Ultrasonics 
Seon) 6 BLeECtronics 


The sound velocity was measured with the pulse echo 
method. A high frequency pulse of 10 MHz and 150 Vip amplitude 
is applied to a piezoelectric quartz crystal of 12 mm diameter 
for 2 usec with a pulse repetition frequency of 100 Hz. The 
transducers for longitudinal sound are x-cuts, for transverse 
I use AC-cut plates. The transducers resonate mechanically 
at that frequency atid send a sound pulse through the helium 
which is reflected, travels back and is picked up again by the 
Same transducer which converts it into an electrical signal. 
The time delay between the initial pulse and the echos is 
measured with the calibrated delayed sweep of an oscilloscope. 
The total path length is obtained from a flight time measure- 
ment in liquid helium and the known sound velocity of liquid 
helium (Vignos and H. Fairbank 1966). I have also determined 
coos velocity at .4,13,K in the pressure range from 20. to 50 
bar by measuring the cell dimensions at room temperature with 
a micrometer and get values which are higher by about 0.2%, 
which is well within the combined errors. 

Figure 19 is a schematic diagram of the circuitry. 


All components except the matching network are commercially 
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FIGURE 19 


Electronics for sound velocity measurement. 
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Cc SOUND TRANSDUCER 


whe) 


available." AS: pulsed oscillator I used part of .@ Sperry 
Ultrasonic Comparator (Sperry Products Inc., Danburry, Conn.). 
The pulse amplifier of that unit failed repeatedly and I 
replaced it therefore with an Arenberg wideband amplifier 
(Arenberg Ultrasonic Laboratory, Jamaica Plain, Mass.) with 

a bandwidth from 0.7 to 70° MHz. This corresponds to a rise 
time of about 0.02 usec which is much smaller than the rise 
time of the pulse generator (0.1 usec). After amplification 
the pulses are rectified and displayed on a Tektronix Model 
556 oscilloscope which is triggered by the Sperry unit. Its 
delayed sweep was calibrated with a crystal controlled counter 
(Transistor Specialities Inc. ,“ Plainview, L.i*) -of- arated 
accuracy OLS parts. 1n 107 and found to agree with the dial to 
Within 0.1%. This calibration was” repeated” after one year 

and no change in accuracy was observed. 

To increase the accuracy of the time measurement I 
read the time from the trigger to the leading edge of as many 
echos as possible (always more than two) and made a linear 
Jeast square fit through a,plot of time vs. echo number.. The 
slope then gives the averaged time between two echos while 
the intercept represents the sum of the delays in the elec- 
tronic circuits probably associated with finite trigger rise 
times. A eee although smaller effect was also seen by 
Vignos and H. Fairbank (1966) who attribute it to a delay in 


the conversion of mechanical into electric signal. It is of 
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the order of 1 or 2 usec in our case, depending on the time 
base setting of the oscilloscope and can therefore not be 
neglected. 

The accuracy of the time measurement is limited to 
Cake see by “the-risetime. of the pulses ‘aridietomovie! by ‘the 
resolution of the dial of the calibrated sweep. Since the 
velocity data used to calculate v are accurate to 0.3% I 
estimate the absolute error in v to be about 0.5%. Using the 
same equipment, the resolution of the time measurement can 
be increased by displaying with high time resolution (0.05 
usec/com) a portion of an unrectified echo with fixed delay 
and observing the phase shift of the echo. The jitter in the 
delayed sweep of 0.01% sets the limit for the sensitivity. 
I found it convenient to record the shape and amplitude of 
the echos during crystal growth. This was done with a boxcar 
integrator (Princeton Applied Research, Princeton, N.J.) and 


a chart recorder (figure 20). 


peace Matching 


The impedance of the quartz transducer at resonance 
is of the order of 1 MOhm in parallel with the electrode 


capacitance C of 10 pF, while the impedances of amplifier and 
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oscillator are around 50 Ohms. The capacitance represents 

a short circuit to the high frequency pulse and a much higher 
oscillator amplitude is required than if C were absent. The 
usual procedure is to put an inductance in parallel and tune 
the resulting LC-circuit to resonate at the mechanical reso- 
nance of the transducer, but this has the disadvantage to 
require critical tuning, which may. not survive during cooling 
of the apparatus and more seriously, it works only for one 
frequency. A solution which solves both this problem and the 
mismatch between oscillator and transducer is found from 
transmission line theory (Skillig (1951), Petco! Liew Ls 
theory says that a transmission line can be used as an im- 
pedance transformer if its length is comparable to the wave- 
length of the electromagnetic wave that it transmits. At 10 
MHz the wave length is about 20 m, and since the coax cable 
between equipment and transducer is 3 m that is the case. The 


complex input impedance ar of a cable with cable impedance 4, 


and length a for a load Z4 and a wave vector k is given as 
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Various graphical ways exist to evaluate this ex- 


pression (Smith (1944)). The first thing to do is to calculate 
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FIGURE 20 


Sound echos in liquid and solid helium during crystal growth 


The elapsed time since the start of growth is given in minutes 
beside each chart recorder trace. Between 0 and 90 min the 
echos correspond to propagation in liquid helium. After 10 
min. a small echo from solid helium is observed and grows 

to reach a maximum amplitude at 240 min. In order to observe 
the small echos in solid helium the amplifier gain was set 
very high. The smaller sound absorption in liquid helium 
then yields echos so large that the amplifier saturates and 
they appear with a clipped off top and constant amplitude on 
the recorder trace. The irregular shape of the multiple 
echos is not well understood but may be due to resonance and 


interference effects in the amplifier. 
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the input impedance of the coaxcable on top of the cryostat 
for the desired frequency from its known length and charac- 
teristic impedance (50 Ohm) and from the known load. The rest 
is then a room temperature problem, namely matching this 
calculated impedance to the amplifier and oscillator. 

This again can be done with transmission lines by 
connecting various impedances along the line from generator 
to load and in parallel to the load and is referred to in the 
literature as double stub matching. Again the place of. 
connection and magnitude of the necessary matching impedances 
can be calculated in principle. In practice the parameters 
of that calculation, namely the generator and load impedance 
and the characteristic impedance of the cable inside the 
cryostat (# 39 copper wire in a stainless steel tube of 0.8 mm 
ID) are not known too well. I have therefore adapted the basic 
layout of a double stub matching network and terminated the 
stub ends with a variable capacitor or inductor and find by 
trial and error the setting of these terminations that yield 


maximum signal into the amplifier. 
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4, RESULTS 


4.1 Sound velocity in hep He 4 as a function of crystal 


orientation. 


The velocity of longitudinal sound for different 
orientations was measured for pressures of 36, 58 and 128 
bar. Table 2 gives total path length 2x1, birefringence 
net TN. (Vos et al. 1967 c), molar volume V and melting 


temperature Tn (Grilly and Mills, (1959), (1962) ) for 


these experiments. 


TABLE 2 
p 2x1 "Ne MO. y Th 
(bar) (cm) (xq Orr) (cm~ /mole) K 


Pout e. 4.5752 0.501°2.77 £0.08 20. 3220.01-1594920-005 
Poms ol. 4.091 U.0L 2085 2.0508) LO 2ero. Oe 2 oe Ol 005 


Meo OSG vie 17 2 0202) 3735 +: 0.08 1723320701 “3992440005 


The first experiments at 58 and 128 bar were done with the 
preliminary bomb. Only two points were measured at 128 bar, 
both giving the same orientation y = 46° and velocity ViiF 
748m/sec. At 58 bar a total of 60 crystals were grown 


(figure 21). The large scatter in the data, which is outside 
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the estimated error in velocity and orientation is probably 
due to an unsatisfactory optical arrangement that could not 
detect polycrystals. The light beam has to pass: 3 mirrors 
between laser and helium crystal and it is likely that it is 
cut off at the edges and does not cover the whole cross. section 
sampled by the sound beam. The criterion of total extinction. 
to establish whether a single or polycrystal is present then 
loses its power and the data suggest that it fails frequently 
in this arrangement. If another criterion is chosen for the 
quality of the crystais namely the number of echos that could 
be detected, the picture improves somewhat. (Black dots in 
Figure 21 are samples that showed more than 6 echos). In any 
case itis clear that the two points at 128 bar have to be 
regarded with caution. 

I designed therefore another sample container with 
a copper bottom that should provide a better nucleation point 
than the quartz window of the preliminary cell and had a 
horizontal optical path which reduced the number of mirrors 
to 2, and thereby increased the success rate for single 
crystals from 5% to 50%. A further improvement would be to 
have windows in the side of the cryostat and eliminate the 
mirrors altogether. This is also the set-up chosen by Lee's 


group (Lee et al. (1969)), and it allows direct visual in- 
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FIGURE 21 


Longitudinal sound velocity in hcp tee V = 19.28 bm? /mole. 


Full circles are from crystals with more than 6 echos. 
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FIGURE 22 


Longitudinal sound velocity in hcp Read % = 20:32 cm?/mole. 
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ELGURE: 23 


Anisotropy of sound velocity in hcp Hees 


All velocities are multiplied with a normalization factor 


to give a longitudinal velocity of 600 m/sec at y = 0. 


NW Nosanow and Werthamer (1965) 
GKW Gillis, Koehler and Werthamer (1968) 


LCH Lee, Crepeau and Heybey (1969) 
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spection of the crystals. 


With this new container I have grown 42 crystals 
(figure 22) at a pressure of 36 bar. For 25 of these I 
determined the orientation. Of the remaining 17 crystals, 10 
showed only one sound echo of an intensity that was smaller 
by a factor 10 or more compared with the other samples. Three 
samples had small spurious, echos with a delay of about 80% of 
the delay for the first echo, which I interpret as a reflection 
from a grain boundary. Since there is a discontinuity in the 
wall of the container at the corresponding distance (the 
window hole pierces through the wall) it may be that crystal 
growth is disturbed at this point. Three samples were lost 
through melting when I tried to anneal them very close to the 
melting temperature. Three more samples were lost for other 
reasons, such as running out of liquid helium in the dewar. 

When it turned out that some of the remaining 25 points 
showed scatter larger than expected from an error analysis, 
their growth record was re-examined and I found that one had 
drifted 1% in velocity during the orientation measurement. 
Three others had the lowest number of echos of all the 25 
samples, namely 3 and I rejected them so that finally 21 points 
remain. These criteria are arbitrary and are only justified 


because the scatter gets smaller as the demands on crystal 
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quality are raised. 

The measured anisotropy of. the longitudinal sound 
velocity can be compared with Lees measurements at a smaller 
density (Lee et al., 1969) (Figure 24) and with two theoretical 
calculations (Nosanow and Werthamer (1965), Gillis et al. 
(1968)). As a guide to the eye, I have drawn a line in 
Figure 22 which represents the calculations of Gillis et al. 
for a molar volume of 19.2 cm?/mole, but multiplied all their 
velocities with 0.915 to get a best fit to the measured points. 
Figure 23 includes in addition the results of Lee et al. and 
of Nosanow and Werthamer, all multiplied by a scale factor to 
give a velocity of 600m/sec for y = 0. I have thereby assumed 
that the anisotropy is the same for different molar volumes. 
Within the limited accuracy of this measurement there is 
quantitative agreement between all four results, with the 
measurements being slightly in favour of the calculations of 
Gillis et al. Also shown in figure 23 are transverse veloc- 
ities, multiplied by the same scale factors. Our results, 
calculated from the elastic constants obtained in section 4.3 
are between the experimental. results of Lee et al. (1969) and 
the theoretical results of Gillis et al. (1968) and have 
qualitatively the same dependence on the angle y. For clarity, 


the data of Nosanow and Werthamer are shown for the directions 
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[0001], [0110] and [0111] SxndYRRy24 0°. 90° ana 46°41°). 

I have measured the transverse velocity and orientation 
in one sample at 36.0 bar and get 319 m/sec at y = 62°. With 
the scale factor chosen for a best fit to the longitudinal 
velocity, this point agrees well with the theoretical and 
scaled down velocity of Gillis et al. (1968) and also with 
the velocity calculated from our elastic constants. 

To my knowledge this is the first observation of 
transverse sound ina He? crystal that was not prepared from 
the superfluid. Growing crystals from superfluid helium is 
thought to be the easiest method and to give the best single 
crystals and from my experience it is more difficult to 
observe transverse sound than longitudinal. This may be the 
reason for the lack of observations of transverse sound above 
29 bar. Lipshultz and Lee (1965) have measured transverse 
velocities in unoriented samples around 25 atm and obtain 
values between 230 m/sec and 315m/sec and Lee et al. (1969) 
finds 225m/sec for the smallest and 309m/sec for the biggest 
velocity in oriented crystals at the same pressure. 

To answer the question whether the observed anisotropy 
an cp fie® is in any way unusual, I have calculated the veloc- 
ity surface of all 16 hexagonal crystals for which the elastic 


constants are available (Landolt-Bérnstein, (1966)). It turns 
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out that there is some correlation between the ratio of. the 
crystallographic axes c/a (Crystal Data, 1963) and the shape 
of the velocity surface. This ratio is obtained from x-ray 
Piemeutron ditfraction and-=is 1.63 for hep we nearly identical 
Bo, the.ratio for close, packed spheres, c/a.= (83)? mo 1663299), 
In crystals that have the same c/a as helium, such as Co 
erro )ymg (1.62) or Ice (1.63), the anisotropy.1s quali-— 
tatively the same with the highest velocity occuring in the 
direction of the c-axis and a minimum around 55°. This is 
also the case for a simple theoretical model of a hexagonal 
solid with harmonic springs between nearest neighbours, for 


which c = 29:11:8:32:8 (Leibfried, (1955)). 


lee Ie 184 38° ° 44 
In crystals for which.c/a deviates much from 1.63, for 
Peemoie in Cd (1.88), Zn (1.86), Apatite (0.73) sor Beryl (1.00), 

the anisotropy is quite different (figure 25) and it appears 
therefore that it is given to a large extent by the shape of 
the elementary cell and-not so much by the exact nature of the 
interparticle potential. Helium seems to be no exception to 
enat rule. 

I conclude this section with a comment on growth in 
meeterred directions. If-the single crystals would grow or 


recristallize from polycrystals in a random way, the number of 


measurements would be proportional to sin y, that is, we would 
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FIGURE 25 


Longitudinal sound velocity in hexagonal crystals. 


All velocities are multiplied with a normalization factor 


to give a longitudinal velocity at y = 0 of 600m/sec. 
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see much more crystals with y = 90° than with y = 0°. It is 
obvious from figure 22 that this is not the case. Growth in 
preferred directions might be induced by irregularities at the 
nucleation point or by the temperature gradient during growth 
and annealing which in our case is perpendicular to the direc- 
tion of observation. If-the c-axis would be preferentially 
along the gradient we again would measure much more crystals 
with y = 90°, which we AG not. All we can therefore say at 
present is that we observe nonrandom distribution of orien- 


tation and that the c-axis is not likely to coincide with the 


temperature: gradient. 


a,2. Sound velocity as a function of molar volume 


The compressibility of solid helium is about 10° times 
higher than that of sodium or about 10> times that of diamond. 
It is therefore possible to reduce the molar volume of helium 
appreciably with modest pressure and to measure the volume 
dependence of the sound velocity or of other properties of 
solid helium. The former is perhaps the most stringent ex- 
perimental test that can be applied to lattice dynamics in the 
long wave length limit. It should be remembered that its 
primary aim is to calculate the energy of the solid, and that 


the volume derivative of the sound velocity involves the third 
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derivative of energy with respect to volume (the first gives 
the pressure, the second the compressibility). 


In figure 25 I compare the measured sound velocity 


along the c-axis (y = 0) and at y = 45° from the c-axis with 
values observed by Lee et al. (1969 b) at a molar volume of 
20.97 com>, Also shown are sound velocities calculated with 
equation 2.38 from the phonon spectrum of hcp single crystals. 
This was obtained from inelastic neutron scattering by 
Peenewrc2 Ct al. (1968). at V.= 21.1 cm? /mole and by Brun et 
my 01969) at V-= 16.0 cm?/mole. The black dots represent 
measurements by Vignos and H. Fairbank (1966) on polycrystals 
or possibly single crystals of unknown orientation. 

The different measurements agree within the error 
limits although the sound velocities of Lee et al. (1969) seem 
to be somewhat low. I have therefore calculated the Debye 
temperature from their elastic constants and obtain eae 
2,58 K which is also lower than Edwards and Pandorfs value of 
26.1 K at that volume, The Bulk modulus KES obtained from 
equation (2.36) is also 10% lower than measured by Jarvis et 
al. (1968). The data of Vignos and H. Fairbank show a smaller 
anisotropy than our measurements. The explanation may be that 
the largest velocities occur in directions around the c-axis, 
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we assume crystal growth with random orientation. In addition 
the lower velocities are associated with large deviations 
between wave vector and beam direction and could therefore 

not have been observed by Vignos (see section 2.3). 

To discuss the volume dependence of the sound velocity 
and to relate it to the volume dependence of the Debye temper- 
ature, the compressibility and the second sound velocity I 
will introduce the Griineisen constant y and begin with a review 


Of Q2ES Various definitions: 


1) Mode Grtineisen constant 


dlog w (q;,s) 
ia;s. a d log V 


This iS a microscopic quantity where w (q,s) is the frequency 
of a phonon with wave vector q and polarization s. V is the 
molar volume. In general Ya,s will depend on orientation, 
polarization and frequency of the phonon, and the molar 
volume of the solid and the temperature. The various macro- 
scopic definitions of the Gruneisen constant will be averages 
in which different frequencies, polarizations or orientations 
will have a different weight, depending on the macroscopic 
quantity. 


Grlineisen's assumption was that in a lowest approxi- 
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mation all modes would have the same volume dependence, in 
which case the various definitions of y will all give the 
same value. In a real solid, this need not be so, nor will 


the Gruneisen constants be independent of volume or temperature. 
2) Gruneisen constant from sound velocity 


In the long wave length limit, where there is no 


dispersion and 


w(q,s) = 9... Vv, 452 
we find 

. — _ Glog. 1 

Ys (ol Kale Pegi 3 ee 


where Me 1s the sound velocity in’direction q.- iv depends 


on direction and polarization s (= 1,t,,t»)- 


3) Gruneisen constant from Debye temperature 


Bets 2109 .8 4 . §4 34 
Yp dlog V 
9, is the Debye temperature at absolute zero. This is the 
quantity usually referred to without specification in the 


Solid helium literature. Expressing 9. in terms of sound 
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velocities 


and differentiating with respect to log V we get 


u =3 -3 
eee NeekVg G@/Eof v_" ao 4.6 


Under the assumption of elastic isotropy and for a typical 


case of v, = : Vi (vy and v 


: being the longitudinal and 


i 


transverse velocities) one gets, 
a (l6y, iz Y,) fRLIS= We 4.7 


Since the transverse modes contribute much more to the Debye 
temperature than the longitudinal, the Debye Gruneisen 


constant will also reflect their dominance. 


4) Grtineisen constant from compressibility: 


dlog v 
Hy li Lae 4.8 
k dlog Vv. 3 } 


where Vy is a quantity with the dimension of a velocity 
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@ is the density, K the compressibility at absolute zero. 

K can be expressed in terms of the elastic constants which 

in turn are obtained from the sound velocities. If we know 
their volume derivatives. Y, Can be calculated exactly. For 


the isotropic case 


ia SP Oa Steg iar a ets, 
ee ei Fey SV, YD 4.10 


The above relations all concern the behaviour of long wave 
length phonons at 0 K and do not depend on any assumptions 


about the equation of state, while the next definition does. 


5) Gruneisen constant from thermal expansion 
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It can be shown by thermodynamic arguments (Edwards and 
Pandorf, (1965)) that nee is independent of temperature and 
equal to Yp if the specific heat .can be written in the form 
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6) Griineisen constant from second sound 


dlog v 
< heelindalt pea 4.12 
i dlog V 3 
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sound velocity (equation 2.39) in such a way that the trans- 
verse velocities are weighted even heavier than in the Debye 
temperature, we would expect that Yrr would be very close to 
Yer but nothing definite can be said before the anisotropic 
case is calculated. 

The experimental information on Gruneisen parameters 
is collected in table 3. We note that both the compressibility 
and the Debye temperature give values for Yp and Ye that are 
lower. at smaller molar volumes. This is in agreement with 
the volume dependence of the sound velocity which also decreases 
with smaller molar volume as can be seen from figure 26, and 
also with the neutron diffraction data at 16.03 om Waoikes which 


shows the decrease of the slope dlog v/ dlog V. 


The Griined sen parameter for longitudinal sound that 
we observe is the same for the two directions of propagation 
that were investigated and is ‘Greer 3.0 2.0.1 for a*vorune or 
20 com>. This is larger than the Debye Gruneisen parameter 


he 2.6 at this density. We can ask how large the transverse 


Vy would have to be to give the correct Debye Gruneisen constant 


Yp and I find, by trying different values of Ye in equation 
4.6 that ea 2.54. I have assumed that y, is the same for 


both transverse branches and independent of orientation. Let 
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TABLE 3 


Gruneisen constants of hcp He’, 


Parameter Density Range Reference 
Yiongitudinal 
(y¥=_0°) #40. 1s 1 99s). Sot ome This work 
Ytransverse 
(y = 45°) jasc: ORM LON re This work 
YD Page wd SO Us Eaten Baa Te repeat E Edwards and 
Pandorf (1965) 
Yp 2.4 Bed eer 3 Dugdale and 
Franck (1964) 
YD Peas Tee 2) Heltemes and 
Swenson (1962) 
Yp Tele FeO OB Seo Leste ee Oe Ahlers (1967) 
Y 24 On - f Base 1 6a 5e8n21 Edwards and 
“ Pandorf (1965) 
Y 2.4 d 2rara wel 6 Dugdale and 
K Franck (1964) 
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us now calculate the volume dependence of the compressibility, 
based on ‘ca 3.0 and Vee 2.594. To.do this i’ calculate the 
sound velocities at two volumes V and V + dV, (using the 
anisotropy of Gillis model), from these the elastic constants 
and the compressibilities from equation 2.36 and finally Vy 
from equation 4.9 for both volumes. I obtain Te Se Which 
compares well with values obtained from measured compress- 
abilities at 20 om? (table 3). The observed volume dependence 
of the longitudinal sound velocity is therefore consistent 
with calorimetric data. 

The work of Jarvis et al. (1968) has shown that the 
assumption that solid helium obeys a reduced equation of state, 
that is C(T,V) = £(T/O (V)) is not justified and therefore 
no simple comparison can be made with the parameter Xa 
derived from thermal expansion measurements. It is also not 
clear why the parameter from the second sound velocity, Yrr = 
2.9, differs so much from Wags but more detailed work on second 
sound, both experimental and theoretical, is needed to discuss 
Baus point.further. 

We finally compare the volume dependence of the longi- 
tudinal sound velocity with theoretical calculations of Nosanow 


and Werthamer (1965) (NW) and of Gillis, Koehler and Werthamer 


(1968) GKW). For simplicity, figure 27 shows only the 
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FIGURE 27 


Theoretical and experimental longitudinal sound velocity 


in theft c-direction. in hep He* (double logarithmic) 


NW Nosanow and Werthamer (1965) 


GKW Gillis, Koehler and Werthamer (1968) 


The experimental curve is the same as in figure 26. 
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velocities in the direction of the c-axis. We see that both 
theories agree well with the observation at a molar volume. of 
20 om?, but also that the volume dependence is not satisfactory 
and leads to discrepancies of. 20% in Vv, at high and low molar 
volume. Stated in other words, the theories predict. Griineisen 
parameters of ‘opae 0.7 (GKW) and sila 1.0 (NW) in disagreement 
with our observed value of Mie 30% 

Which approximation in the theory is. mainly responsible 
for this is not clear at present and different possibilities 
are discussed by Gillis et al. (1968): 

In both calculations 3-particle short range corre- 
lations were neglected. We would expect that at smaller 
volumes, where the motion of the particles is more restricted, 
these would play a more important role. They have been in- 
cluded in a calculation of sound velocities at.different 


: by de Wette et al. (1967) and change 


densities for bcc He 
the Gruneisen parameter Yy from.1.0 (only- 2-body terms in- 
cluded) to 1.3.. They also found. that the cluster expansion 
converges less well for the hcp phase so that the inclusion of 
3-body terms might lead to more significant improvements in 
this phase. 


Nevertheless it is thought that the main cause of the 


discrepancy is the particular analytic form (equation 2.13) of 
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the correlation function £, which has the advantage of pro- 
ducing a quickly converging cluster expansion. A better 
correlation function would improve the theory at small volumes 
and at the same time lead to a better convergence, which is 
needed for the hcp phase. Work in progress on this problem 
has been reported at the Aspen Quantum Crystal Conference 
(1967) by Horner and Guyer and Nosanow. 

Another approximation to be checked is the inter- 
particle potential which may not be well represented by a 
Lennard-Jones potential (Munn and Smith, (1965)). However 
Since the large zero point motion produces an average over 
this potential its detailed behaviour is probably of minor 


importance. 


23) The: elastic constants-of solid helium 


Of the five elastic constants. of a hexagonal crystal, 
only. four (Crys Cy3r C337 Cay) determine the longitudinal 
sound velocity; Cy> Can only be determined from the pure 
transverse mode. Of these four constants only two can be 
obtained directly from the present set of measurements with 
sufficient accuracy: C), from the longitudinal velocity in 


the basal plane and C33 from the longitudinal velocity in the 


c-direction. 
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To determine the remaining constants we may use the 
tongitudinal velocity: at.y = pepe published data for the 
compressibility and a relation between four of the elastic 
constants that holds for all hexagonal crystals provided they 
do not change their shape under hydrostatic pressure. 

Solid helium satisfies this condition very well at 
least up to a pressure of 1100 bar and probably beyond as can 
be seen from the observed variation of the ratio c/a of the 
crystallographic axes with pressure (table 4). The most 
accurate determination of c/a comes from birefringence 
measurements because these are sensitive to deviations from 
the ratio for close packed spheres, c/a = 1.63299, while 
neutron and x-ray diffraction experiments measure c anda 
separately. The c/a ratios calculated by Vos et al. (1967 b, 
1967 c) using the theory of Kronig and Sonnen (1958) change 
my. bY ~@ parts-in 107 for a volume change from 12.5 to 20.6 
cm?/mole. This implies that the linear compressibility (the 
relative decrease in the length of a line when the crystal is 
subject to hydrostatic pressure) must be the same in the c- 
direction and in the a-direction in solid helium. The ex- 


pression for the linear compressibility in a hexagonal crystal 


is (Nye (1964), pel146) 
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where y is the angle between c-axis and the direction of the 


line. The condition 8(0) = 8 (90°) leads to 


Cc. +c - cCc - Cc = 0 4.14 


and the volume compressibility K becomes now 


K_ t=o3/ (C3, + 2c 


3 13? 


From published values for K and the previously determined 
constants Ci and C33 we obtain C13 from equation 4.15 and 
then C15 from equation 4.14. The three constants Cia [13 and 
C33 together with the longitudinal velocity at y = 45° finally 
yield Cay from Peer on Jie to's 

Values for the isothermal compressibility in the 
density range of interest have been published by Edwards and 
Pandorf (1965) and by Jarvis et al. (1968). Since the temper- 
ature dependence of Edwards and Pandorf's isothermal com- 
pressibilities is in conflict with the measurements of Jarvis 
et al. and also with our measurements on the temperature de- 
pendence of the sound velocity, we consider Jarvis data more 


reliable. We find by interpolation of. their measurements to 
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the densities of our measurements for the compressibility at 
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0K, K, = 0.309 x Oe” om? /ayn ati Vy =:20732 cm?/mole and 


: om? /dyn AR gn ge Net | cm?/mole. The corresponding 


O2 218 x10: 
velocities at y = 45° are 517 m/sec and 596 m/sec, res- 
spectively. A correction would have to be applied to all 
velocities to get the elastic constants at 0 K because of the 
variation of the velocities with temperature. This correction 
of about 0.2% however is small compared with the errors in 

Vi and is therefore ignored. The full set of elastic constants 
obtained in this way is shown as a function of molar volume 

in figure 28. Table 5 gives numerical values together with 
constants directly determined by Lee et al. (1969) from their 
measurements of longitudinal and transverse velocities. Also 
shown are elastic constants obtained with equation 2.33 from 
published theoretical sound velocities. We estimate the 
accuracy of our constants to be + 4%. 

These constants can now be used to calculate the 
velocities of the transverse modes and the velocity of second 
sound from equation 2.33 and 2.39 (Table 6). The one measure- 
ment of the transverse velocity at 20.32 cm?/mole is 2.0% 
lower than the value of 327 m/sec calculated from the elastic 
constants, which we consider satisfactory agreement. 


Interpolating second sound velocities reported by 
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FIGURE 28 


Theoretical and experimental elastic constants of hcp He* 


vs. molar volume. 


NW Nosanow and Werthamer (1965) 
GKW S6illis et al.‘ (1968) 
LCH ‘Lee et al. (1969) 
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TABLE 6 
Longitudinal, transverse and second sound 


Velocity for different orientations y 


psec? 0:30 iem~ /nidle 
0° 600 m/sec 267 m/sec 267 m/sec 157 m/sec 

10 595 Z67 275 159 

20 579 266 296 163 

30 554 264 322 166 

40 528 263 345 167 

50 510 261 352 167 

60 508 259 Bo E65 

70 Si / 258 305 162 

80 525 25./ 278 L577 

90 528 256 267 154 

Vee=2 19.28 om’ /mole 

Y ab Vel Vi2 be ie 
aa i ll Scale ee RO Rc Seren A at 

on 687 m/sec 307 m/sec 307 m/sec 180 m/sec 

10 681 306 316 es 

20 663 305 338 187 

30 637 302 Ges 190 

40 608 299 Sho RE 190 

50 588 295 398 188 

60 585 292 380 186 

70 593 289 348 182 

80 602 Zot 319 her ay 


90 ai eee 287 307 174 
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Ackerman and Guyer (1968), we find 130 m/sec at 20.32 cm?/mole 
and 170 m/sec at 19.28 emmole. It is not clear at present 
why the measured velocities are smaller than the values shown 
in table 6, but it should be remembered that equation 2.39 
assumed elastic isotropy which is certainly not the case in 
solid helium. Qualitative agreement exists between the 
angular variation of the second sound velocity observed by 
Ackerman and Guyer (1968) and our data. 

I have also calculated the Debye temperature at 0 K, 
Oo" by numerical integration of equation 2.35 and obtain Cras 
Boece O.6 K at 20.32 com> /mole and om = 33,2 27056) Raat. ao 
cm?/mole in good agreement with calorimetric values of 28.3 K 
and 32.5 K, respectively, reported by Edwards and Pandorf 
12965). 

It is remarkable that one of the Cauchy relations, 
Cag = C37 (equation 2.40) is satisfied very well, especially 
at the higher densities. This is the relation that depends 
on the assumption of central forces. The other, c,, = 3C} 51 
which requires in addition that each atom is a center of 


inversion symmetry, does not hold however, as is expected in 


eehcp ‘lattice’ 
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4.4 Sound velocity as a function of temperature 


Using the high resolution method for measuring small 
velocity changes described in section 3.31 we have measured 
the temperature dependence of the sound velocity at a volume 
OL 20.42 cm? /mole (p = 34.2 bar) between 0.7 K and the melting 
temperature of 1.9 K. We used-a He? cryostat of conventional 
design and a high pressure cell similar to sample chamber 1 
described in section 3.12, except that the quartz window was 
replaced by a polished bottom of high conductivity copper to 
improve crystal growth. Blocking of the high pressure line 
during crystal growth was prevented with a heater wire inside 
the line. The heater current was switched off during the 
velocity measurements, which therefore give values at constant 


volume. The data can be fitted to an equation 


4 4.3 


(Py = 57 eeCl= Claes AWLO) xt Da) AUG 


6T AllOK} 7 in m/sec) 


The high absolute value of the velocity indicates that the 
propogation direction was close to the c-axis and the sample 
was probably a single crystal. 

No theoretical estimates exist at present for the 


variation of the sound velocity with temperature and we there- 
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fore discuss the temperature dependence of the compressibility 
velocity Vv, (T) (equation 4.9) which should give an idea of the 
order of magnitude of the variation of the longitudinal 
velocity. 

The temperature dependence of the isothermal compress- 
ability Kp (T) was found by Jarvis et al. (1968) who measured 
(dp/dT) with a strain gage technique for different pressures and 
obtained Kn at the melting temperature from the work of Grilly 


and Mills (1962). K,(T) can then be obtained from numerical 


integration of the thermodynamic relations 


00 
ey lee 
oT p eke) 7 


where ‘ 4.17 


K,(T) has also been calculated from the specific heat by 
Edwards and Pandorf who assumed that solid helium obeys a 


reduced equation of state. With this assumption they get 


v1 (oy = tp ((y.+1) (U-U_)-y,TC_) 4.18 
Kn (T) - Ky = RY, Yp 0? pry: ° 


U- U. is the temperature dependent part of the internal 
O 


energy and is obtained from the specific heat Co. The 
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FIGURE 29 


Ratio of specific heat at constant volume, C to specific 


she 


heat at constant pressure, .C_, vs. temperature. 


Vv =120.4 cm-/mole. 
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temperature variation of K,, calculated in this way is only 


A's 
about half as large as the more directly obtained value of 
Jarvis et al., who also show that their compressibilities can 
not be represented by a reduced equation of state. We there- 
fore think that equation 4.18 is wrong if applied to hep He* 
Edwards compressibility would also result in an increase of 
velocity with temperature, while we observe a decrease. 

At nonzero temperature we have to distinguish between 
isothermal compressibility K,, (temperature T constant) and 


adiabatic compressibility K, (entropy S constant). The 


S 


relationship between them is 
K./Kn = CH Mp 4.19 


where GE and Sanare heat capacity at constant volume resp. 


pressure. Thermodynamics gives the ratio of the specific heats 


as 
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This equation now includes only quantities that are directly 
accessible to measurement. eC ORE EY > 1 is given in Jarvis (1968) 


paper. For a reduced equation of state we would get 


2 
K oe ae 
: ear ey Ste Garena) A232 


We see that K,. - K, 18 proportional to rt ine this 


T 
approximation (figure 29). 


At a molar volume of 20.4 ea 


K, (T) increases by 
about 3% between T = 0 and the melting temperature Gime hk 
and Kn (T) is larger by about 2.4% than K, (T) at this temper- 
ature. 

Whether the sound velocities are given by the 
adiabatic or isothermal compressibility (or rather: elastic 
constants) depends on the frequency (Bhatia, (1967), p. 268). 
For low frequencies the wave length is sufficiently large that 
the compressed (or heated) and rarefied (or cooled) regions in 
a sound train are far enough apart that no heat conduction 


occurs between them. The compression is therefore adiabatic. 


If the frequency is higher than a characteristic frequency 
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ge noe 4 vi/k 4.23 


(k is the heat conductivity; factors of the order unity are 
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ignored) this is no longer the case and the isothermal constants 
determine the longitudinal velocity. The transverse velocity 
is obtained in any case from the isothermal constants since no 


volume change is involved in their propagation. 


In the intermediate region the velocity Vy 


lo Mee 4.24 


2) = 
Vv = (0.K_) 
k Ss B ea 
(Kp /K,) “r We 


The expression for Spe is equivalent to 


W = T ey. EN 4.25 


where a and X are the relaxation time and mean free path for 
umklapp processes. Values for \} have been calculated by 
Mezhov-Deglin (1965) from his thermal conductivity measurements, 
they depend on temperature as Te in the temperature range 
that we have investigated. In our case the sound velocity was 
measured at a frequency of 5 MHz which is smaller than WO, 

down to about 0.8 K. At this temperature we would therefore 
expect a velocity change of the order of 5 (C{0. 8K) /(C,, (0.8K) 
=) ny io Unfortunately this is just at the limit of 
resolution of our velocity measurement but should be easily 


observable with more sophisticated equipment (Abraham et al. 


(1969)). Associated with the velocity change is a peak in 
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FIGURE 30 


Temperature dependence of the sound velocity in solid helium 


at V = 20.4 eni/moLe (estimated). 
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sound absorption (Guyer (1966)). 

Due to the uncertainty of our knowledge of the relax- 
ation time ae and also because we have neglected factors of 
the order unity in equation 4.23, the transition temperature 
of 0.8 K is quite tentative at present, and in addition nothing 
is known about the temperature dependence of the transverse 
velocities. However, the general behaviour of the longitudinal 
velocity V, as a function of temperature should be quite 
Similar to the one of Mi. shown in figure 30. 

In figure 31 we show our measured velocity MS and also 
Vy calculated from the temperature dependence of the 
adiabatic compressibility obtained from the data of Jarvis et 
al. and from Edwards and Pandorf. The change Av, = Vv, (T) = 
Vv, (0) from Jarvis et al. is approximately proportional to 


ie in the temperature range T = l...2 K, in good agreement 


of the observed dependence of Av, a hoe For a Debye solid, 
where 90(T) is independant of temperature, Av, would be pro- 
portional to 4, Tiehncp He® this approximation applies for 


Peery s0D°. O or T < 0.6 K in our case. 
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FIGURE 31 


Longitudinal sound velocity in hcp He* vs. temperature. 

V = 20.42 om?/mole. The solid lines represent the temperature 
dependence of the velocity calculated from the temperature 
dependence of the adiabatic compressibility based on measure- 


ments. of Edwards and Pandorf (1965) and Jarvis et al. (1969). 
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5. CONCLUSIONS 


It 1S possible to grow large single crystals of 
solid helium at pressures between 30 and 130 bar and to 
measure their orientation using the birefringence technique. 
The longitudinal sound velocity in these crystals is largest 
in the direction of the hexagonal c-axis and has a minimum in 
a direction 55° away from the c-axis. This is typical 
for many other hexagonal crystals. 

Extrapolations of the measurements to higher and 
lower. densities agree with sound velocities obtained from 
neutron scattering experiments and also with unpublished 
sound velocities in crystals grown from superfluid helium 
by Lee et al. (1968). The variation of the longitudinal 
velocity with volume can be expressed by a mode Gruneisen 
parameter * ee 3.0 + 0.1, which is higher than the value 
found from specific heat measurements, RO= 2.6.) THiS: isnot 
unexpected and similar to the situation found in other. 
materials. 

Measurements of the longitudinal velocity in the c- 
direction, in the basal plane and in an intermediate direction 
are sufficient in solid helium to determine all five elastic 
constants if the compressibility is known. The Debye temper- 


atures calculated from these elastic constants were found to 
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agree with the calorimetric Debye temperatures within the 
combined errors of the two measurements. 

For a molar volume 20.42 cm?/mole the longitudinal 
velocity decreases by 0.25% between T = 0.6 K and the melting 
temperature of 1.9 K, which is consistent with the increase 
of the compressibility with temperature, observed in strain 
gage experiments. 

All these results lead to the conclusion that from 
an experimental point of view solid helium is not different 
from any other solid in its elastic behaviour. 

The measurements allow for the first time a detailed 
comparison of the predictions of the theory of quantum crystals 
with experiment. At our densities, good agreement between the 
magnitude and the anisotropy of the longitudinal velocity is 
found, but present theory predicts a much weaker volume. 


dependence than found experimentally. 
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7. APPENDIX: COMPUTER PROGRAMS 


All computations were done. on a IBM/360 time sharing 
system using the language APL/360 (Falkoff and Iverson 


(1968)). 
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SOUND PROPAGATION IN UEXAGONAL CRYSTALS 
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IS THE DEBYE TEMPERATURE AT ABSOLUTE ZERO, 

D SPECIFIES THE NUMBER OF SECTIONS INTO WHICH 
THE UNIT SPHERE IS DIVIDED FOR THE NUMERICAL 
INTEGRATION. 2 IS ACCURATE TO ABOUT 0.03 


PERCENT FOR D+10. 


ett 


, aaa iat ANEIAC ORV POV OW 
AL caItnOASOHA vag a LAT AAS ti 
«2 ) BHOLTOAALA WOUTMD SWABAIG BOD DW 


QWA esxn-9 HERNNR a 


oe oe Oar 


- QTUNTSHOD 
e*«¥D\o 40 BTL M 13 Se 


{ 


(LIMO. ate AOLPOWYA BOF “QanlUORA) sata 

; ye AM ML YTPLOOISY ANWLCUTTOWOS 

3 Vaz 2AM WL YRLOOIEN. SRABVAMAST of 

' oak: ogayt WL LRIDOD AR AURAVLMART 6S 


anser SAUL CUTLASS 40 aanoe ML WOVERIVES 9 > OATH 
TANSOWAV AM MORE WARS (Ai 


r 


MORE MASE ques. BRAANEMANT 19 a ae ATIEG 1 


AOTOEY WOLTASTAAIOS FO 40 SayANsG WH OLEREY NG ray 
MOLTIEALG QOUTAETRALOS TAUT MONT 


ctv Nav Wi Qaiey HOTPOWT peemmnane ett" xy at cca 


ay h Badsd ‘ 


ra AWKOTLOWYA BIOBA ant 
DRO, BAAADAG w Xe) 


Ltd 


C3] 


LF 


VVZCOIV 

Z2*+G VL C 

G+OGt180 
2+((C01]+0C05])x(10G)*2)+(¢CL414+C0L5]1)«(20G) «2 
Z*100x((Z+G PH C)#2xCL6])*0.5 


VVT1LOIV 

Z+G VT1 ¢C 

G+oG?180 
Z+100x(((0.5x(CL1J-cL£2])«(10G)*2)+CL5]x(20G)*2)+cL 
6])*0.5 


VVT2LOIV 

Z2<G VT2 C 

G<0G+180 
Z+((CL13+C05])x(10G)*2)+(CL4]+CL5])x(20G)*2 
Z*100x((2Z-G PH -C)+#2xCL6])*0.5 


VPALOIV 
Z<G PHC 
Z+(((C0£11-C0L5)])*2)«x(10G)*4)+((CL4)J-CL5])*2)x(20G)x* 


4 
Z<+(Z+2x((10G)*2)x((20G)*2)x((CL1)-cl5J])x(CL5)-cl4ul)))+ 
2x((C031+CL5])*2))*0.5 


VDELTALLOIV 

Z2*+G DELTAL C3H3M3N3A3E 

H<+(CL6I*x((G VL C)#100)*2)-cL5] 

M<100G#180 

N<200G4180 

A+(H- (cCL4j- CLSJ)xN*2)+ #MxMxCL31+CL5 J 

F<~208((( (107 304) *4)+Mx2)+(( (20 30A)*4)+ Nx2))* 
G.. 5 

Z<+(Hx30F)+(CL61x((G VL C)#100)*2) 

Z«( 30Z)x180+01 


VDELTATICOIV 

Z2<+G DELTAT1 C3E3H 

E« 20100G#180 

H+(CL61x((G V71 C)+100)*2)-c[5] 
Z<(~30(Hx30F)+#(CL6]*x((G Vf1 C)#100)*2))x180+01 


ca (908) (L2TO4E# IOD4 (Gs (DOE )*CERIDATET 
&.o*(Ld JD*S4 (9: Bo ‘Pr Be: 


’~ 


s(posyx(subteme -forp) ton coor ico CABLED | 


oe 
4% 


x 
2 


Ra 


- 3 S)x(fen- cJo Yet a(S ae 
sala #8 @ 04 (S#(C2 39406 30) Ix 
ais 4 “yu yaaa 
VIDA 

io 2 eH WENO ARTA | 
tan aegynasn ce: NO) dat (a9) 7 


Lar 


\ 
p 
«4 


i 

j 

Te. 
ae 2 


here 7 


~ 2 
2 ba 

yy le” a0) 5 = 
a << se 


a 


. a 
© 


08 
ty yay 
C2Ipele oo rane etece2 88) 
« ( (Set (Pe (NOE: “eoy))4+¢ w (me (noe "e i 
| Ba ae: ty A 
aor ‘en(0or@ 8 dyoete 1994 (OEXB) 


co i 
g * 


138 


VDELTAT2COIV 

GCG DELTAT2 Ci:H3sM:N3;A:3E 

HeCUCL61x((G V72 ¢)+100)*«2)-cL[5] 

M<100G4180 

N<200G#180 

A<«(H-(CL4I-CL5)])xW*2)4+NxMxCL31+C05] 

E«” 202((((10 304) *4)+M*2)+(( (20° 30A)*4)3N*2))* 
0.5 

Z4<«(Hx30H)+(CL6]x((G V22 C)+100)*2) 

Z+(° 30Z)x180401 


VDELLLOIV 

6G DELL C3HsA3M3H 

H<«(CL61]*x((G VEZ C)4100)*2)-CL5] 
M+100Gi180 

N<«200G4i180 

A+ 30(H-(CL41]-C05])xN*2)+MxNxCL314+C05] 
Z<«( 20((10A)xM)+(20A)xWV)x180401 


VHEXACLOIV 

Z<N HEXA C 
Z+(((904NV)4+1),8)p0 
Z0311]«0,(1(90N))xW 
20324-20357): VE-C 

WL SarreZitsij VTi. ¢ 
Sh l+2(211 VT? ¢ 

vid Oo ae ON 

Z203;51+2031] DELTAL C 
2136]+20£3;1] DELTAT1 C 
203;71]«2Z031] DELTAT2 C 
203;81]+20£3;1] DELL .C 
2Z(1,((90#N)+1)3; 5 6 7 8]+0 
Z;.171<0 


VDEBALOIV 

Z<+D DEBH C3:G3;M 

Z<0 

M<1i 

G«(~ 454D)+Mx904D 
Teta (GS (| (UB C) MR VET O)AG Vela C))*-3)x100G+180 
+((M<M+1)<D)/3 

Z<0. 36273433x( 2xDxCL6]+0ZxC[ 7] )*#3 


BEL | 2 ty. 
: ‘ FV j 
<oym e 
a Behe ¥i Fy 
cancteetttiogs 
& 


‘ cap teeealeneceed pote) 
#((sete (etnies “os )) 4 Ses (#4 (ROE™ un iaae 


Ae 
Se(OOLe(D STV >))xfaIn#c3 : dod. 
“ee vires S087 eR 


sets ene 


we 
ae Wats Asked. 4130 9 +e: | 
f2J3- (sx(oore( av D))xE aJ2 oe on 
| gatenoes aa ‘(ted 
a Teate ibe eneataciete a Oe eh ai 
a bene COR te A CN NA bad 


sm o 


“Daca, centunme)? 
wx ( (0400) 1). 0+ 


_ 


aoe | 

2. stataad Loe] 

> 130 0 

one #8 2 p(te(Meoe EIS 
9° qi r¢ Pee’ fn} er biat. Ss & 


4 / ee 
- A 7 ee iy —_- 7h 
m Se Pa Sy “groe: MM ae : ; 
a rsdoorx(e-#((9 cay ae 8 ae 


] | CAAT +.) ieee © bY 


EXAMPLE 


GeSeo 2:94 1.31 7.1 1.403 0.1970 4.003 


20075..0 


DT i ge SI 


woe VT! C 


10 DELTAL 


10 DELTAT 


10 DELTAT 


10 DEEL C 


20 DEBH C 


5 HEXA ¢ 


600.3 
098. 9 
594.7 
36759 
578.6 
567.3 
554.4 
941,0 
928.0 
oop i OF 4) 
509.7 
506.9 
508.14 
oid. 9 
346.7 
921.3 
925.1 
0 oe) 
928.4 


C 


i 


2G 


266.9 
266.8 
266.6 
266.2 
265.7 
265, 0 
264.3 
20 iso 
2.6 25.0 
261, 7 
260.8 
Phe 
209060 
258.3 
257.6 
25st 
256.7 
25645 
256.4 


266.9 
268.9 
eo hee 
20049 
Boo 8 
308.7 
322.3 
335.1 
345.4 
hae Oy | 
352.4 
346.6 
335.4 
320.7 
304.9 
289 we 
ie 37 
205 57 
266.9 


° ° 


° 


° 


° 
OFOFN OMA GG vj 


° 


ay 
ONWOAMAWOMDE 
° ° 


af 1S Se ee eet el eae Foe Eom ae See ae a ce 
DOOFPPRPNNNNNNFPFPEH OOO 
e ° 


° ° ° s es ° ° 


° 


Wee -Oe oe fo, es Se 


° °o e ° ° ° ce ° 


fs fs bo fs 


SOP WOFFEFENOFDWDONNFPUOATUNO 
Ss 
ONPNAIIN AYOYINHNOPPRPRaAUHE OO 


° 


°o ° ° ° ° ° ° s 


° 


33 


oret.o eon. te tea 


OPMrOAtrh ot oO 


ge » 


‘ 
i 


vn Ss ae 
“er? a t Hae? s 


nt eee 


abpecose ae 
a4 


, Sis bates; 
ae met 
Mau pee bia tet 32 
io 7 a 
- STATISG he ~~ 
jheeer 
‘-, 2 6fes Gee “a 


ae bs Fe 
hil eect i 


Gipsy ne aan 


ce a 8° BeORe 50S 
; z a 
rr Bey tog “2 
ie ll pe aailhes ia - 
A os hub ms aah ; 7 ~ 
= *) HM 


w, ‘ye. 
we c ag 


B29970 


